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ABSTRACT 

We have developed a one-dimensional photochemical and thermochemical kinetics and diffusion 
model to study the effects of disequilibrium chemistry on the atmospheric composition of "hot Jupiter" 
exoplanets. Here we investigate the coupled chemistry of neutral carbon, hydrogen, oxygen, and ni- 
trogen species on HD 189733b and HD 209458b, and we compare the model results with existing 
transit and eclipse observations. We find that the vertical profiles of molecular constituents are 
significantly affected by transport-induced quenching and photochemistry, particularly on cooler HD 
189733b; however, the warmer stratospheric temperatures on HD 209458b help maintain thermochem- 
ical equilibrium and reduce the effects of disequilibrium chemistry. For both planets, the methane and 
ammonia mole fractions are found to be enhanced over their equilibrium values at pressures of a few 
bar to less than a mbar due to transport-induced quenching, but CH4 and NH3 are photochemically 
removed at higher altitudes. Disequilibrium chemistry also enhances atomic species, unsaturated hy- 
drocarbons (particularly C2H2), some nitriles (particularly HCN), and radicals like OH, CH3, and 
NH2. In contrast, CO, H2O, N2, and CO2 more closely follow their equilibrium profiles, except at 
pressures < 1 microbar, where CO, H2O, and N2 are photochemically destroyed and CO2 is produced 
before its eventual high-altitude destruction. The enhanced abundances of CH4, NH3, and HCN are 
expected to affect the spectral signatures and thermal profiles of HD 189733b and other relatively cool, 
transiting exoplanets. We examine the sensitivity of our results to the assumed temperature structure 
and eddy diffusion coefficients and discuss further observational consequences of these models. 
Subject headings: planetary systems — planets and satellites: atmospheres — planets and satellites: 

composition — planets and satellites: individual (HD 189733b, HD 209458b) — 

stars: individual (HD 189733, HD 209458) 



1. INTRODUCTION 

The detection and characterization of extrasolar 
planets is one of the most exciting and fastest growing 
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fields in astronomy ( e.g.. j Cha rbonneau et al.' '200"^ 
.Marley et a l. 2007; Ud ry et al.l 12007: Dcming & Scmm 
I2OO9I: iBaraffe et al.l 12010 1 Of the handful of tools 
available for characterization of exoplanet atmospheres, 
transit and eclipse observations have proven to be 
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invaluable for infe rring properties of both th e planet and 
its host star (e.g., 
Henry et al.|[200a 



Charbonneau et all [2000 . [2002. 2005]; 
Seager fc Sasselov'12000^ 'Brown .200!l: 
Brown et all 120011: Burrows et al. 2001; Hubbar d et al l 
2001 1 : iDeming et al. 2005b; Ki iutson et all I2007L 
200a iMadhusudhan fc Seag er. ^QQ^" The wave- 
length dependence of the signal during transit and 
eclipse observations can allow identification of atmo- 
spheric constituents on "h o t Jupiter" exoplanet s (e.g. , 
Charbonneau et al. '2002'; 'Vidal-Madiar e t all 120031: 
Tinett i et al. 2007; GriUmair ct al. 2008; Swain et al.l 
2008bL l2009at IMadhusudhan fc Seageii l2009| ). allowing 



a better understanding of the chemical and physical 
processes operating in these exotic atmospheres. 

Hot Jupiters are often divided into two classes, 
those that appear to have stratospheri c thermal in- 
versions, and those that do not (e.g.. iFortnev et aD 



2008at iHubenv et al.l 
Burrows et al.l 



iHarrington et al.V 



.200.1 ■ ■ 

■ .20071 120081: iKnutson et al.l 

Spiegel et aLl 120091: IMadhusudhan fc SeagelT 

the ultimate reasons for these dif- 

still under investigation (see also 

Spiegel et al . 2009: Zahnlc ct al. 2009; 



Although 
ferences are 



20Ta rYoiTdin fc Mitchelll 120101: iFressin et al.l 



Machalek et al. I 120081 120101: IMadhusudhan fc Seaged 
2010al) . exoplanets with thermal inversions appear to 



have some source of high-altitude absorption, whereas 
exoplanets without thermal inversions are apparently 
missing this stratospheric absorber, pe rhaps because 
it is tied up in condensed phases (e.g., IHubenv et al.l 
120031: IFortnev et al.ll2008al: iBurrows et al.ll2008D or is de- 
stroyed by photolysis and/or other mec hanisms related 
to ch romospherically active host stars (jKnutson et al.l 
l2010h. The observations discussed in the the above 
citations suggest that HD 209458b has a strong strato- 
spheric thermal inversion, whereas HD 189733b does 
not. These two transiting planets, which have been 
the most extensively studied and characterized to date 
because of their bright parent stars and large transit 
depths, become convenient endmembers for our study of 
the chemistry of the two main classes of hot Jupiters. 

Constraints on the composition and thermal struc- 
ture of HD 189733b and HD 209458b have been 
provided through (1) the identification of Na, H, 
and other atomic neutral and ionic species in 
the p lanets' upper at mospheres ( Charbonneau ct aLl 



200l iVidal-MadiaTe ral. 2003, 2004, 2008: Ben-Jaffe 
2007, 2008; Ballcstcr ct al. 2007; Snellen ct al. 200" 
Redfield et al. 2008 : .Sing et al... 2008a; .Ehrenreich et all 
20081: iLan gland-S hula et al.l 120091: iLinskv et al.l 120101 : 
Ben-Jafrelll2010t iLecavelier des Etangs et al.l 120101 ) , ~(2) 



the identification of molecular species such as H2O, CH4, 
CO2, and CO in the ~1 bar to 0.1 mbar atmos pheric 
region (Tingtti et al.l 120071; I Fortnev fc Marleyl 120071: 
Burrows et al. 2007, 2008: Grill mair et al.ll2008l: iBarman 
2008; Swa in et al.ii2008aiib[ .,2009a"'bl'2010^%eaulie u et al. 
2008, 20091: iDesert et al .1 12009: Madhusudhan fc Seager 
2009), (3) observational and theoretical inferences re- 
garding the vertical and horizontal temperature struc- 
ture, including the suggested presence of a stratospheric 
temperature inversion on HD 209458b but not HD 
189733b, and the identification of l ongitudinal tempera- 
ture v ariations and phase lags (e.g., Knutson et al. '2007|, 
120091: iCooper fc Showman. .2006; .Fortncv ct al.. .2006ai. 



2008al [2010| iRichardson et al l 120071; IBurrows et al l 



2007 



2009; 



2008F 



Cowan et al.l l2007t _ Showm an et al.l 120081 
200a .2010aj) . and (4) 
the inferred presence of hazes, clouds, and Rayleigh 
scatterers th at can affect atmospheri c spectral be- 
havio r (e.g., Charbonneau et al.l 120021 : IFortnev et al.l 
2001 IFortnev! 120051: llro'eT al. 2005; P ont et all 12008 : 



Lecavelier des Etangs et al Il2008a >b,: .Sing et al.ll2008allb 
20091: iShabram et al.l \20m . Note, however, that the 



data reduction and observational analyses are often diffi- 
cult due to low signal-to-noise, potential variable stellar 
activity of the host star (including star spots), and/or 
poorly understood systematic instrument effects (e.g., 
iGibson et al.ll2010l) ; spectral models are also poorly con- 
strained, making analysis of the observations difficult. 
Various groups do not always agree, ev en when ex- 
amin ing the same or similar da tasets (cf. [De sert et al.l 



20091 vs. [Beaulieu et al.l 120081 vs 



Swain ct al.ll2008bl v s. Sing et al.ll2009l vs. IGibson et~aLl 



Ehrenreich et al.£^200 



2010t iStevenson et al.ii201Q vs. lBeaulieu et al.ll201lD . An 



additional problem with analyses is the lack of reliable 
molecular line parameters for hot bands of gases other 
than H2O, CO, and CO2. Even if the data reduction 
and analysis were not inherently difhcult, interpreta- 
tion of thermal infrared data can be hampered by the 
fact that molecular bands can manifest as either absorp- 
tion or emission features, depending on the atmospheric 
thermal structure, creating a degeneracy between tem- 
perature profiles and compositions (e.g., 'Burrows et al 
12008; Madhusudhan fc Seager 20 01120103; Tinetti et al 
I2OIO; Fort nev et al.l 1201 0: O'Do novan et al.ll2010f ). De- 
spite these complications, transit and eclipse observa- 
tions represent our best current means for characterizing 
the atmospheres of hot Jupiters. Theoretical models are 
needed to interpret these observations. 

Existing chemical models for transiting exoplan- 
ets tend to fall into two groups: thermochemical- 
equilibrium models and photochemical models. 

Thermochemical equilibrium is assum e d for most 

mode ls (e.g., Bvutows fc Shard 119991: iSeager et aTl 
'2000!; ILodde rs fc Fcglcv 2002; S harp fc Burrowsl 
1 20071: IBurrows et al. 2007. ,2008: Fortnev et al.l I2008al: 
iVisscher et al.l l2010al) . although potential departures 
from equilibrium are sometimes explored through 
evaluations of simple time-scale argument s (e.g 
Fodders fc Feglev 2002; Fortnev et al. 2006a, 2008 



Visschcr ct al. 2006, 2010a; MMhusudhan_fc_ Scage: 
20111 ) or through more cor nplex models with p arame- 
terized chemical kinetics (jCooper fc ShowmanI l2006t) . 
Thermochemical equilibrium is a reasonable starting 
point for exoplanet composition predictions, but the 
strong ultraviolet flux incident on close-in transiting 
planets like HD 189733b and HD 209458b can drive 
the observable upper regions of the atmosphere out 
of equilibrium. Moreover, some of the existing trans- 
port and kinetics time-scale argu ments have known 
flaws (for f urther discussion, see lYung et al.l 11988 
Smith 1998; Bczard ct al. 2002"; 'Vissch er et aI\\20TOh 



Moses et al. 2010; Visschcr fc Moses 20 111) and are often 
misused, leading to questionable predictions of quenched 
disequilibrium abundances. 

The second group of chemical models exploits knowl- 
edge of Jovian photochemistry to predict the conse- 
quences of the large ultraviolet flux incident on ex- 
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oplan et atmospheres (iLiang et aLl 120031 



latter models, 



iKoskinen et al 



_ _ _ , _ 2007. 2009; 

Zahnle et al.i i200a~l2011l: pLine et alj 12010 ). Of these 



2004 iGarc ia Mun og 120071: [Koskine n et al 



YeMj2(M\ IGarda Mufioa (pOOl . and 
• (120071 . 12009h focus on the uppermost re- 
gions of the atmosphere, i.e., the thermosphere, iono- 
sphere, and exosphere and do not co nsider what goes on 



deepe r th an 1 microbar. In contrast. iLiang et alj (|2003l 
120041 ) and lLine et al.l ()2010[ ) explore the photochemistry 
of the deeper stratosphere and troposphere using thermo- 
chemical equilibrium and/or transport-quenched disequi- 
librium as a lower boundary condition for their models; 
however, their focus on low-temperature reactions (for 
the former two studies) and their lack of fully reversed 
reactions (all three studies) prevents them from extend- 
ing the kinetic models into the thermochemical regime at 
temperatures > 1000 K, and their assumptions can result 
in artificial constant-with-altit ude profiles fo r the the ma- 
jor species in the troposphere. iZahnle et al.l ()2009l [20Tll) 
have published the only models to date that have at- 
tempted to bridge the thermochemistry and photochem- 
istry regimes with a single model; however, their reliance 
on isothermal atmospheric models can result in some 
unrealistic vertical profiles and species abundances that 
will not be present with more realistic atmospheric struc- 
tures, and it is not clear whether their reaction scheme 
will reproduce equilibrium at high temperatures. 

We therefore have developed a unique one-dimensional 
(1-D) thermochemical and photochemical kinetics and 
transport model to circumvent some of the problems 
described above. The model has the ability to tran- 
sition seamlessly from the thermochemical equilibrium 
regime in the deep troposphere of giant planets to the 
"transport-quenched" disequilibrium and photochemical 
regimes in the upper troposphere and stratosphere to 
illuminate how chemical processes in each regime com- 
bine to influence the observational properties of giant- 
planet atmospheres. In an initial set of studies, we ap- 
plied the model to the deep troposphere of Jupiter to 
ensure that the model satisfies observational constraints 
on a well-studied planet within our ow n solar system 
(|Visscher et al.ll2010bl : iMoses et all 120101 ). We now ap- 
ply the model to the atmospheres of HD 189733b and 
HD 209458b, as representatives of the two classes of hot 
Jupiters — those with and without stratospheric thermal 
inversions. Our main goal is to more realistically pre- 
dict the vertical profiles of the major observable carbon-, 
oxygen-, and nitrogen-bearing species on these close-in 
transiting planets to aid the analysis and interpretation 
of primary transit and secondary eclipse observations. 
By including both thermochemical and photochemical ki- 
netics in a traditional 1-D chemistry /diffusion model, we 
can avoid some of the shortcomings and pitfalls of previ- 
ous chemical models — we do not need to correctly iden- 
tify in advance the mechanism and rate-limiting steps for 
transport quenching of species like CO, CH4, HCN, and 
NH3, we can avoid other uncertainties in the time-scale 
arguments such as the correct length scale to use for the 
transport time scale (e.g.. Smith 1998), we do not need to 
restrict the model to any particular pressure region, and 
we can consider realistic atmospheric structures rather 
than isothermal atmospheres. 

As we will show in the following sections, disequi- 
librium processes like transport-induced quenching and 



photochemistry can perturb the composition of HD 
189733b and HD 209458b away from equilibrium. The 
effects of these disequilibrium processes are more sig- 
nificant for cooler exoplanets like HD 189733b than for 
warmer exoplanets like HD 209458b, but there are po- 
tential observational consequences for both planets. In 
particular, transport-induced quenching and photochem- 
istry can enhance the abundance of atmospheric con- 
stituents like CH4, HCN, NH3, and C2H2, and the pres- 
ence of these disequilibrium molecules will affect spec- 
tral and photometric behavior. We compare the model 
results with existing transit and secondary eclipse ob- 
servations and discuss in detail the major disequilibrium 
mechanisms affecting the abundance of carbon, nitrogen, 
and oxygen species on HD 189733b and HD 209458b. 

2. MODEL 

We use KINETICS, th e Caltech/JPL 1 -D photo- 
chemistry/diffusio n code (lAllen et al.l [T98lt see also 
lYun g et al.l [1981 [G ladstone et al.l 119961: IMoses et al.l 
2000allbl. 120051 : ILiang ct al. 2003, 2004; L ine et al.ll20'Tnfl 

to calculate the vertical distributions of tropospheric and 
stratospheric constituents on extrasolar giant planets. 
This fully implicit, finite-difference code solves the cou- 
pled mass-continuity equations as a function of pressure 
for each species: 



'dt 



dz 



P, - L,, 



(1) 



where is the number density (cm~'^), $i is the flux 
(cm^'^ s~^), Pi is the chemical production rate (cm~^ 
s~^), and Li is the chemical loss rate (cm"'^ s~^) of 
species i, all of which are explicit functions of time t 
and altitude z. The flux transport terms include both 
molecular and eddy diffusion, the latter parameterized 
by an eddy diffusion coefficient K^z- The calculations 
are allowed to evolve until steady state is achieved, and 
a convergence criterion of one part in a thousand is used. 
We consider 198 atmospheric levels with a vertical reso- 
lution ranging from 30 grid points per atmospheric scale 
height in the deep atmosphere to 3 grid points per scale 
height in the upper atmosphere. A finer grid is needed 
in the deep atmosphere to accurately track the quench 
points for some of the disequilibrium species. The model 
extends from the deep tropospheres of HD 209458b (bot- 
tom boundary at 2754 K, 1 kbar) and HD 189733b (2600 
K, 1.67 kbar) to well past the homopause region into 
the thermosphere (top boundary at pressures less than 
10^^° bar, exospheric temperature assumed to be 12,000 
K). The upper boundary is chosen at a high-enough al- 
titude such that we encompass the photochemical pro- 
duction and destruction regions for the main molecular 
species H2, H2O, CO, CH4, CO2, N2, NH3, and HCN for 
plausible assumptions about K^z, and such that assump- 
tions about upper boundary conditions do not affect the 
results. Although the models extend into the thermo- 
sphere, our intent is not to specifically model thermo- 
spheric chemistry. Instead, the thermospheric tempera- 
ture profile and escape boundary condition are include d 

^[^ia([2ool . 



in an ad hoc manner based on the results of 



IGarcia Munod (|2007[ ). and IKoskinen et al 
below for further details) 



(|2010ir 



see 



We use free-convection and mixing-length theories 
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Fig. 1. — Eddy diffusion coefficient profiles adopted for our nom- 
inal models, as estimated from rms v ertical velocity profiles ob- 
tained from the GCMs of IShowman et al.. (2P09j) (see text). For 
both planets, the black line represents the profile for the dayside- 
average model, and the gray line represents the profile for the 
terminator-average model. The high-pressure value of Kzz = 10^" 
cm'^ s~^ is an estimate of the eddy diffusion coefficient in th e adia- 
batic region from free-convection and mixing-length theory llStonel 
119761) . and K^z's are assumed to be constant at altitudes above 
the top level of the GCMs. 



(e.g.. lStonelll976HFlasar fc Gieraschlll977[) to define the 
eddy diffusion coefficients K^z in the deep, adiabatic por- 
tion of the troposphere and estimate K^z profiles through 
the bulk of the rest of the atmosphere from the root- 
mean-square (rms) vertical velocities derived from global 
horizontal averages at a given pressure l e vel fr om the 3- 
D GCM simulations of IShowman et all (|2009[ ). i.e., by 
assuming Kzz = w{z)L{z), where i(;(z) is the horizon- 
tally averaged global rms vertical velocity from the GCM 
simulations and L[z) is assumed to be the atmospheric 
pressure sca le height H{ z) (but could be some fraction 
of -ff(-z), see lSmithlll998D . Note that this procedure pro- 
vides only a crude estimate of Kzz, and it would have 
been preferable to have calculated Kzz from the eddy 
vertical velocity times the eddy displacement, but this 
information was not available from the GCMs. In partic- 
ular, the procedure may overestimate Kzz in the ~10-200 
bar radiative region, where the vertical motion tends to 
consist of small-scale wave oscillations. In the future, it 
would be useful to include passive tracers into the GCMs 
to directly and rigorously calculate the rate of vertical 
mixing. Figure [T] shows the eddy diffusion coefficient 
profiles adopted for our nominal models. However, given 
the above caveats, we consider Kzz profiles to be essen- 
tially free parameters in the models, and we examine the 
sensitivity of our results to reasonable variations in the 
adopted Kzz pro files. The mo l ecular diffusion coefficients 
are described in iMoses et al.l (|2Q00al) . 

Temperature profiles are not calculated self- 
consistently within the KINETICS code but are 
instead required as inputs to the model. The tempera- 
tures adopted in our model are shown in Fig. [21 Most 
of these profiles d e rive from the 3-D GCM simulations 
of IShowman et ahl (|2009i) for the region from 17 bar to 
a few /zbar or the 1-D models of iFortnev et al.l H2006al 
I2010D for the region from 1000 to 1.0 x 10"^ bar. 
Because the nitrogen species in particular quench at 



very high temperatures, we needed to extend the GCM 
profiles to deeper levels [i.e., to re ach temperatures of at 
least 2500 K). The 1-D models of IFortnev et"all (|2006al 
120101 ) already include the deep adiabatic region, but the 
GCM pro files do not. For HD 189733b, the 1-D Mtt" 
profiles of IFortnev et al.l ()2006al 2010) that assume the 
incident stellar energy is redistributed over the entire 
planet compare very well with the GCM theoretical 
thermal profiles in the 1-100-bar region, so we extended 
the GCM profiles to deeper pressures by simply con- 
necting smoothly to the 1-D profiles. For HD 209458b, 
the 1-D "47r" profiles are somewhat warmer than the 
GCM profiles in the 1-100-bar region, and we shifted 
the 1-D profiles by -60 K before connecting the GCM 
profiles to this deeper adiabatic region. At low pressures 
(high altitudes), we extended the profiles by assuming 
either an isothermal atmosphere or a high-temperature 
therin ospher e, using the th e rmosp heric profiles of lYelld 
()2004D and iGarcfa MunHg ()200<1) as a guide for the 
latter. Once the pressure-temperature profiles were 
established, the rest of the background atmospheric grid 
(e.g., densities, altitudes) was derived through solving 
the hydrostatic equilibrium equation. 

HD 189 733b orbits its host star at a distance of ^0.031 
AU (e.g., ISouthworthI l2010f) . The cool but chromo- 
spherically active HD 189733 has an effective temper- 
ature in the ra nge -5050-5090 K (jBouchv et al.l 120051: 
IKnutson et al.l l2010) and has typically been classified 
in th e literature as a Kl V-to-K2 V main-sequence star 
(e.g., iBarnes et ai1l2010l : iShkolnik et al...2008). although 
the H ipparcos catalog lists it as a G 5 V (jPerrvman et al.l 
119971: see also iMontes et all 1200 It) . In any case, HD 
189733 is less luminous than the Sun, and HD 189733b 
is not one of the hottest of the transiting exoplanets that 
have been discovered to date, despite its small orbital 
semimajor axis. HD 189733 is part of a binary star sys- 
tem, but the M-dwarf companion to HD 189733 is ^^200 
AU away and should not significantly affect the flux inci- 
dent on HD 189733b. Because well-calibrated ultraviolet 
spectra for HD 189733 are available for only a specific 
range of wavelengths (e.g., iLecavelier des Etangs et ahl 
we set up a normalized 1-AU spectrum using the 
Hubble Space Telescope (HST) STIS spectra of epsilon 
Eridani (a K2 V star) from the CoolCAT database (jAversI 
I2005D for the 1150-2830 A region, after correcting for 
stellar distance. At the upper end of this wavelength 
range, the fiux from epsilon Eridani is much less than 
for typical G stars like our Sun at the same stellar dis- 
tance, whereas the fiux at the shorter EUV wavelength 
end of this range compares well to main-sequence G stars. 
We therefore use the 1-AU solar fiux divided by 10 for 
wavelengths longer than 2830 A, and we use the 1-AU so- 
lar flux for low-to-average solar conditions as an analog 
for our 1-AU-normalized HD 189733 flux for wavelengths 
below 1150 A. The photochemical calculations are then 
performed assuming that HD 189733b is located 0.031 
AU from this K2V-analog star. The host star for HD 
209458b is a GOV stellar type (Hipparcos, Ferryman et 
al. 1997), and we simply use the solar flux as an analog 
for HD 290458, and assu me HD 209458b o rbits at a dis- 
tance of 0.047 AU (e.g., ISouthworthI IMol) . The stellar 
zenith angle 6 is fixed at 48° for the dayside-average mod- 
els — where (cos 6) —2/3 {6 48°) is the projected-area- 
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Fig. 2. — Tempe r ature profiles for our models of HD 189733b (left) and HD 20945 8b (right), based on tf i e gen eral circulation models 
of IShowman et al] II2009I ) (solid lines) and on the 1-D radiative-convective models of IFortnev et all l l2006al . 120101) (dotted lines), for an 
assumed solar-composition atmosphere. For HD 189733b, the models include a projected-area-weighted profile averaged over the dayside 
hemisphere (to best represent the planet's apparent disk during the secondary eclipse), profiles averaged over the eastern and western 
terminators (to best represent the ingress and egress regions during the transits), a terminator-average profile, and 1-D profiles calculated 
assuming the incident stellar energy is distributed over the entire planet (curve labeled "4pi") or over the dayside hemisphere only (curve 
labeled "2pi"). For HD 209458b, the models include a straight (i.e., not projected-area) average over the dayside hemisphere, along with 
eastern and western terminator averages, a terminator average, and an average over the nightside hem isph ere. Most pro files have been 
given an ad hoc high-temperature thermosphere at high altitudes based on the models of I Yelld II2004I ) and IGarcia Muncj^ (I2007l ) . Also 
shown are pressure-temperature curves (dashed gray lines) at which the major nitrogen species, N2 and NH3, and the major carbon 
species, CO and CH4, have equal abundances in thermochemical equilibrium. To the left-hand side of these CH4 = CO and NH3 = N2 
curves, the reduced species (CH4 and NH3) dominate, while the more oxidized species (CO and N2) dominate to the right. Note that 
both planets are within the CO and N2-dominated regimes over much of the atmosphere for all temperature profiles except for the coolest 
HD 189733b western-terminator-average profile. The gray dot-dashed lines represent the approximate condensation curves for the major 
silicates enstatite (MgSiOa) and forsterite (Mg2Si04), based on iVisscher et al.i i2010ai) . A color version of this figure is available in the 
online journal. 

weighted average of the cosine of the stehar zenith angle 
over the planetary disk at secondary eclipse conditions 
— and B is fixed at ~84°-86° for the terminator models. 
Multiple Rayleigh scattering by H2 and He is considered 
in the model using a Feau trier radiative-transfer method 
(jMichelangeh et ailflOOl . 

We consider the neutral chemistry of 90 carbon-, 
hydrogen-, oxygen-, and nitrogen-bearing species: the 
model contains H, H2, He, C, CH, excited singlet ^CH2, 
ground-state triplet ^CHs, CH3, CH4, C2, C2H, C2H2, 
C2H3, C2H4, C2H5, C2H6, C3H2, C3H3, CH3C2H, 
CH2CCH2, C3H5, C3H6, C3H7, C3H8, C4H, C4H2, 
C4H3, C4H4, C4H5, C4H6, C4H8, C4H9, C4H10, C5H3, 
C5H4, C6H2, C6H3, C6H4, C6H5, non-cyclic CeHg, 
CgHe (benzene), O, 0(iD), O2, OH, H2O, CO, CO2, 
HCO, H2CO, CH2OH, CH3O, CH3OH, HCCO, H2CCO, 
C2H3O, CH3CHO, C2H4OH, HO2, H2O2, N, N2, NH, 
NH2, NH3, NNH, N2H2, singlet biradical H2NN, N2H3, 
N2H4, CN, HCN, H2CN, CH2NH, CH3NH, CH2NH2, 
CH3NH2, CH2CN, CH3CN, C3N, HC3N, C2H2CN, 
C2H3CN, NO, NO2, N2O, HNO, HNO2, NCO, HNCO, 
and CH3NO. Metals, rock-forming elements, and sulfur 
and phosphoru s species are not co nsidered, nor is ion 
chemistry. See IZahnle et al.l (j2009D for a discussion of 
the effects of sulfur photochemistry on exoplanet com- 
position. The reaction list derives originally from the 
Jupite r and Satu rn models of [Gladstone et al. (1996), 
iMosesI (fl996[) . and lMoses et"all |i995a .b., .2000a. h. ,2005i) . 
although extensive updates that account for high- 
temperature kinetics have been included based largely 
on combustion-c hemistry literat ure ( e.g., Baulch et all 
[T99l[T99l [20051 : 1 Atkinson et al.l[T997ll20od : ISmith et aLl 



[2000'; 'Tsans' 198?, T99ll IDean fc Boz^ [20001 ). See 

iVisschcr et al. (2010b) for further discu ssion of the hy - 
drocarbon and oxygen kinetics and Mos es et al.l (|2010f ) 
for further discussion of the nitrogen kinetics. The model 
contains ^^1600 reactions, with the rate coefficients of 
^800 of the reactions being taken from literature values, 
and the remaining '-^^800 reactions being the reverse of 
these "forward" reactions, with the rate coefficients of 
the reverse reactions being calculated internally at each 
pressure-temperature point along the grid using the ther- 
modynamic principle of microscopic reversibility (e.g., 
kfor/krev = K^q, whcrc kfor IS the rate coefficient for the 
forward reaction, k^ev is the rate coefficient for the re- 
verse reactio n, and Kf^g is the equilibr ium constant of the 
reaction; see iVisscher fc Mosesl(20Hl for further details). 
The rate-coefficient expressions for the forward reactions 
(and their literature sources) are provided in Table SI of 
the Supplementary Material. The rate coefiicients for the 
reverse reactions for the pressure-temperature conditions 
of certain of our exoplanet models can be found in the full 
model outputs als o presented i n the S upplementary Ma- 
terial. Note that iMiller et all (|2OO90 demonstrate that 
detailed balance assuming microscopic reversibility is ex- 
pected to be quite accurate, even for complicated chem- 
ical reactions that involve multiple, interconnected po- 
tential wells. 

We a ssume a 1 x protosolar composition at mo- 
sphere ([Loddersl [200l [20091; iLodders et all [20091 ) in 
thermochemical equilibrium for our initial conditions 
for all the con stituents, using the NASA CEA code 
([Gordon fc McB ridc 1994) for the equilibrium calcu- 
lations. Thermodynamic parameters are taken from 
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Gurvich etall (I1989-1994D the JANAF tables (IChasd 
19981) ■ iBurcat fc Ruscid (|2005D . or other hterature 
sources. Note that although we do not consider rock- 
forming elements and their resulting effects on the chem- 
istry of oxygen, we do assume that ^20% of the oxygen 
has been sequ estered along with silicates and metals (see 
lLoddersl[200l in the deep atmospheres of HD 189733b 
and HD 209458b, thus depleting the available oxygen 
above the silicate clouds. As such, our "solar" abundance 
of oxygen is slightly less than that used by many other 
groups, and the species profiles themselves should not be 
considered accurate within and below the silicate cloud 
condensation regions due to our neglect of the effects of 
rock- forming elements (see Fig. [2]). 

For boundary conditions, we assume zero fiux at the 
top and bottom boundaries so that no mass enters or 
leaves the system. The lower boundary in our mod- 
els is at a high-enough temperature that thermochem- 
ical equilibrium will prevail, and a zero-flux boundary 
condition is reasonable. We also checked cases with a 
fixed-mixing-ratio lower boundary condition for all con- 
stituents, set at their equilibrium values, rather than a 
zero-fiux lower boundary condition and found no differ- 
ence in the results. The upper boundary in our model 
is at a pressure of a few times 10~^ /ibar, which typi- 
cally corresponds to a planetary radius less than 3 times 
the 1-bar radius Rp in exoplanet hydrodynamic models 
(and in so me cases much less than 3i?T» d e pending on the 
model, see ISniHOOl iTian et al.l [20051 Ictarcia Munozl 
|2007| ). At these high altitudes, the hydrodynamic mod- 
els and model-data comparisons show the atmosphere of 
HD 209458b to be escaping, with an overall mass loss 
rate that is low enough that exoplanets like HD 209458b 
would lose less tha n ^1% of their mass ov er their lifetime 
(e.g..[Lecavelier des Etanes e t al. 2004; T ian et al.ll2005l : 
IMurrav-Clav et al.l I2OO90 As such, we also test cases 
with an upper escape boundary condition for species like 
atomic H, with a typical upward flux of ~ (6— 60) x 10^^ 
cm-2 s-i (e^ 10^-10^ g s"! for HD 209458b at a radius 
of 3R „, see r^lle 2004. Garcia M ufloz 2007, Ti an et al.l 
[20051 and IKoskinen et al.l I2010l) . or we investigate a 
range of possible upper boundary escape fluxes from 10^*^ 
to 10^^ cm^^ s^^, but we flnd that this upper bound- 
ary condition has no effect on the stratospheric or lower- 
thermospheric results. However, we have not included a 
hydrodynamic wind connecting to this escaping region, 
which is the main reason the escape flux has little ef- 
fect on our model stratosphere. Hydrodynamic m odels 
(I Yelld 120041: fTian et al.l[2005[: IGarcia Mimo3[2007l) show 
that the atmosphere should be close to hydrostatic equi- 
librium below 2-3i?p, with the base of the planetary wind 
usually bein g located near the 1-n bar level (see also the 
discussion m iKoskinen et~alll2010[ ). Thus, escape from a 
hydrodynamic or thermal wind is not expected to sub- 
stantially influence the composition of the stratosphere, 
and hydrostatic equilibrium is a reasonable assumption 
in our models. However, if an atmospheric bulk wind ac- 
tually dominates down to the base of the thermosphere 
(or deeper) for any particular exoplanet, then hydrody- 
namic flow and atmospheric escape could conceivably af- 
fect the stratospheric composition on such exoplanets. 

3. RESULTS 



Our model results suggest that the three main chemi- 
cal processes — thermochemical equilibrium, transport- 
induced quenching, and photochemistry — all oper- 
ate effectively in the atmospheres of HD 189733b and 
HD 209458b. Thermochemical equilibrium dominates 
at pressures greater than a few bars, transport-induced 
quenching can dominate for some species in the ~1-10~'^ 
bar region, and photochemistry can dominate at pres- 
sures less than ~10~^ bar, except when the stratosphere 
is very hot. All three processes combine to influence the 
vertical profiles of the important observable constituents 
on HD 189733b and HD 209458b. As such, all three 
processes should be considered in investigations of the 
chemistry and composition of extrasolar giant planets. 
Below we describe how the different disequilibrium pro- 
cesses affect the results, and we discuss the dominant 
mechanisms that control the abundances of neutral car- 
bon, oxygen, and nitrogen species in our thermochemical 
and photochemical kinetics and transport models. The 
column abundances for some of the interesting species in 
our models are shown in Table [T] and full model outputs 
are supplied in the Supplementary Material. 

3.1. Equilibrium vs. transport-induced quenching 
vs. photochemistry 

Fig. O illustrates how the individual disequilibrium pro- 
cesses affect the mole fractions of some species in our 
dayside-average HD 189733b and HD 209458b models. 
The dashed curves show the thermochemical-equilibrium 
solutions. Note that CH4 dominates over CO, H2O dom- 
inates over CO, and NH3 dominates over N2 at equilib- 
rium in the deepest regions of both atmospheres, but 
the dominance reverses at higher altitudes such that CO 
and N2 are the dominant carbon- and nitrogen-bearing 
species in equilibrium at pressures less than ^10 bar (i.e., 
over most of the observable portion of the atmosphere; 
see also Fig. [5]). At the highest altitudes {i.e., in the 
thermosphere), atomic species C, N, and O become the 
dominant neutral carbon, nitrogen, and oxygen carriers. 
Although H2O remains an important carrier of oxygen 
at all pressures on these planets, thermochemical equi- 
librium models predict that the mixing ratios of CH4 and 
NH3 will drop off dramatically with increasing altitude, 
such that CO > CH4 and N2 > NH3 in terms of the 
column abundance at 1 bar and above. Thermochemical 
equilibrium predicts an even smaller column abundance 
for species like HCN and C2H2, such that these species 
would be at most minor constituents on HD 189733b and 
HD 209458b if their atmospheres were in equilibrium. 

The situation changes significantly when vertical 
transport is considered. The dotted curves in Fig. [3] 
show the results of models in which thermochemical 
kinetics and vertical transport are operating but in 
which the planet receives no photolyzing ultraviolet 
radiation from the host star, i.e., so that we can examine 
the infiuence of transport-induced quenching separately 
from that of photochemistry. At deep pressure levels in 
these atmospheres, temperatures are high enough that 
energy barriers to the kinetic reactions can be overcome; 
reactions proceed in both the forward and reverse 
directions, and equilibrium is maintained. Our kinetics 
models, which have fully reversed reactions, reproduce 
equilibrium predictions in these high-temperature, 
high-pressure regions. However, as a gas parcel is lifted 
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TABLE 1 

Column abundances above 1 bar 





nominal K^^ 


nominal K^^ 


K,, = 10'' 


K,, = 10" 


nominal Hzz 


nominal Kzz 




terminator 


dayside 


dayside 


dayside 


dayside 


terminator 


Species 


HD 189733b 


HD 189733b 


HD 189733b 


HD 189733b 


HD 209458b 


HD 209458b 



H 


2.1 


X 


1020 


9.3 


X 


1020 


1.2 


X 


1021 


7.9 


X 


1020 


H2 


1.0 


X 


1026 


1.0 


X 


1026 


1.0 


X 


1026 


1.0 


X 


1026 


He 


1.9 


X 


1025 


1.9 


X 


1025 


1.9 


X 


1025 


2.0 


X 


1025 


C 


5.6 


X 


1014 


4.1 


X 


1015 


1.2 


X 


1015 


6.1 


X 


1014 


CH3 


2.0 


X 


lOi'' 


2.3 


X 


1017 


7.2 


X 


IOI6 


3.7 


X 


IOI8 


CH4 


5.5 


X 


1021 


1.5 


X 


1021 


4.4 


X 


1020 


2.4 


X 


1022 


C2H2 


6.6 


X 


1016 


3.0 


X 


1016 


2.8 


X 


1015 


4.6 


X 


10i« 


C2H4 


1.8 


X 


1016 


6.3 


X 


1015 


5.6 


X 


1014 


1.6 


X 


10i« 


C2H6 


2.4 


X 


lOii 


5.4 


X 


1013 


4.9 


X 


IOI2 


1.3 


X 


1016 


C3H2 


9.4 


X 


1014 


9.2 


X 


1014 


7.1 


X 


1014 


2.5 


X 


1015 


C4H2 


4.3 


X 


1012 


1.5 


X 


lOlO 


1.5 


X 


109 


6.1 


X 


1013 


CeHe 


6.8 


X 


1012 


4.5 


X 


lOii 


8.1 


X 


10" 


2.5 


X 


10" 





2.5 


X 


1014 


7.7 


X 


1014 


8.5 


X 


1013 


5.8 


X 


1014 


OH 


1.6 


X 


1015 


8.3 


X 


1015 


7.9 


X 


1015 


1.3 


X 


IOI6 


H2O 


4.6 


X 


1022 


4.2 


X 


1022 


4.1 


X 


1022 


6.7 


X 


1022 


CO 


5.1 


X 


1022 


5.4 


X 


1022 


5.5 


X 


1022 


2.9 


X 


1022 


CO2 


1.6 


X 


1019 


1.2 


X 


1019 


1.2 


X 


1019 


9.9 


X 


IOI8 


HCO 


1.7 


X 


1014 


4.7 


X 


1014 


4.8 


X 


1014 


2.5 


X 


1014 


H2CO 


1.0 


X 


IOI6 


1.0 


X 


IOI6 


1.1 


X 


IOI6 


5.6 


X 


1015 


CH3OH 


1.8 


X 


1013 


7.3 


X 


lOlO 


8.7 


X 


IOI2 


6.1 


X 


1013 


N 


6.4 


X 


1014 


2.3 


X 


1015 


5.6 


X 


1013 


9.8 


X 


1014 


N2 


7.3 


X 


1021 


7.5 


X 


1021 


8.2 


X 


1021 


4.9 


X 


1021 


NH2 


8.4 


X 


1015 


2.8 


X 


IOI6 


2.5 


X 


1015 


9.1 


X 


IOI6 


NH3 


1.7 


X 


1021 


1.2 


X 


1021 


1.0 


X 


1020 


4.1 


X 


1021 


HON 


1.8 


X 


1020 


1.4 


X 


1020 


5.0 


X 


IOI8 


2.6 


X 


1021 


H2CN 


1.3 


X 


1013 


1.8 


X 


1013 


8.1 


X 


10" 


3.6 


X 


1014 


CH2NH 


1.9 


X 


1015 


9.3 


X 


1014 


3.9 


X 


1013 


1.8 


X 


IOI6 


CH3CN 


6.2 


X 


1014 


1.6 


X 


1014 


1.9 


X 


1012 


4.8 


X 


IOI6 


NO 


1.7 


X 


1014 


4.3 


X 


1014 


8.5 


X 


IOI2 


5.2 


X 


1014 


HNCO 


4.5 


X 


1015 


5.0 


X 


1015 


4.4 


X 


1014 


9.0 


X 


1015 



8.9 


X 


1021 


2.3 


X 


1026 


4.5 


X 


1025 


2.1 


X 


1015 


1.7 


X 


1017 


2.1 


X 


1020 


2.7 


X 


1015 


2.4 


X 


1014 


7.7 


X 


10" 


3.7 


X 


108 


1.9 


X 


105 


1.0 


X 


10-2 


3.5 


X 


1015 


1.5 


X 


1017 


9.6 


X 


1022 


1.3 


X 


1023 


2.3 


X 


1019 


3.4 


X 


1015 


2.3 


X 


IOI6 


9.4 


X 


IOI2 


2.1 


X 


1015 


1.9 


X 


1022 


2.2 


X 


IOI6 


1.8 


X 


1020 


1.3 


X 


1019 


4.1 


X 


IOI2 


6.4 


X 


1013 


1.2 


X 


1012 


4.6 


X 


1014 


1.0 


X 


1015 



2.2 X 1022 

2.3 X 1026 

4.5 X 1025 
7.0 X IO15 
1.9 X 1017 
1.2 X 1020 

3.4 X IO15 

1.6 X 10i4 

3.2 X 10" 
7.9 X 106 

3.3 X 104 

1.2 X 10-4 
X IOI6 
X IO17 
X 1022 

1.3 X 1023 

1.9 X 1019 

5.0 X IO15 
2.3 X IOI6 

7.3 X IO12 

7.1 X IO15 
1.9 X 1022 
3.9 X IOI6 
1.6 X 1020 

1.4 X IO19 

5.5 X IO12 

5.6 X 1013 

7.7 X 10" 
7.6 X IO15 
1.1 X IO15 



6.1 

5.5 
9.4 



Note. — Model column abundances have units of cm" ; X^^'s have units of cm 3- . "Terminator" refers to a terminator- 
averag e and "dayside" refers to a dayside-averagc temperature profile for the planet in question from the GCMs of lShowman et al.l 
j200gf) . Our "nominal" K^z profiles are shown in Fig. m 



to cooler, higher altitudes, energy barriers can begin 
to become insurmountable over typical atmospheric 
residence time scales, and exothermic "forward" re- 
actions can become favored over their endothermic 
reverses. When the transport time scale drops below the 
chemical kinetic conversion time scale between different 
molecular species, the species can become quenched 
such that their mole fractions become "frozen in" at 
values representative of the quench point (i.e., where 
the tw o time constants are equal). iPrinn &: Barsheivl 
(|1977| 1 first developed this concept analytically to 
explain the unexpectedly large CO abundance in 
Jupiter's upper troposphere. Based on CO and other 
observations, transport-induced quenching appears to 
be ubiquitous in substellar objects, operating on our 
solar -syst em giant pla net s (e. g. 



T977t ILewis fc PrinnI [19801: 



1981 



Prinn fc Barshav 



Prinn fc Olaguer 



ILewis fc FeglevI Il984t iFeglev fc PrinnI 1 198 

Feglev fc Lodderd 11991 'Lodders fc FegleV '20021 
Bezard et all ^ l2002t FVisscher & Feglev 2005; 

Visscher et"al1 l2010bt 'Moses et al. 2010^1, on brown 



dwarfs (e.g.. iFeglev fc Lodders 1996: NoUetaLl 
Griflnth fc Yelld 119991 ISaumon et al.l I200a 
20061: ILeggett et all 120071: IHubenv fc B urrowsl . 
Geballe et al.l 120091: lYamamura et al.l 12010: Kin g et alJ 
20101 iVisscher fc MosesI I2011D . and most likely on ex- 



trasolar plane t s, as w el l (e.g.. iCoooer fc Showmanll2006l: 
'Fortnev et al.' '2006a; 'Burrows et all 120081: iLine et al.l 
1010; Madhusudhan fc Scagcr 120111 and references 
therein). Although both CO and CH4 will quench when 
the kinetic interconversion between the two species 
becomes inefficient, the quenching is more obvious for 
the species that is not expected to be abundant in 
equilibrium (i.e., CO on the cooler giant planets and 
brown dwarfs; CH4 on the warmer, highly irradiated, 
close-in transiting giant planets). 

The dotted curves in Fig. [3] show that transport- 
induced quenching can greatly increase the predicted 
abundances of several key species, including CH4, NH3, 
and HCN, over expectations based on thermochemical 
equilibrium (dashed curves) . The quench point for a ma- 
jor species like CH4 is caused by the breakdown of effi- 
cient kinetic interconversion between itself and another 
major species (CO in this case) and results in an ob- 
vious departure from the predicted equilibrium curve. 
Once an abundant species like CH4 becomes quenched, 
its quenched abundance can affect other minor species 
like C2Ha, that might not have reached their own quench 
points yet; these species may continue to react with 
the CH4 at higher altitudes in an equilibrium with the 
quenched CH4 until they too reach their own quench 
points. Thus, minor species may depart from equilibrium 
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Fig. 3. — Mole fractions for several species in our models for HD 189733b (left) and HD 209458b (right) for different assumptions about the 
chemistry and transport. All HD 189733b models assume a projected-area-weighted dayside average thermal structure with an isothermal 
extension at high altitudes (i.e., with a cold thermosphere), and all HD 209458b models assume a dayside-average thermal structure 
with a hot thermosphere (see Fig. [2] for the temperature profiles for these models). The dashed lines are for a purely thermochemical 
equilibrium model (no kinetics or transport), the dotted lines are for a model that has thermochemical kinetics and transport but no 
incident ultraviolet photons (i.e., no photochemistry), and the solid lines represent our full model with thermochemical and photochemical 
kinetics and transport all operating. The eddy diffusion coefficient Kzz is fixed at a constant value of lO' cm^ s~^ in the transport models. 
A color version of this figure is available in the online journal. 
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Fig. [3] ( continued) . 
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Fig- \S\ ( continued) . 

when a major parent molecule quenches, but their mole 
fractions may not become constant until they themselves 
quench at higher, cooler altitudes. As such, the mole 
fractions from the full kinetics and transport model can 
have complicated behavior that is not easily predicted 
from simple time-constant arguments using a single rate- 
limiting kinetic reaction. Hydrogen cyanide (HCN), for 
example, maintains an equilibrium with H2 and with 
the quenched abundances of CH4 and NH3 until its own 
quench point when the interconversion of NH3-N2-CH4 
shuts down, and the resulting mole-fraction profile is 
complicated. In addition, the quenching may not be 
abrupt, which can affect the mole-fraction profiles and 
make abundances more difficult to predict analytically. 
Ammonia, for example, has two quench points in the HD 
189733b model. Ammonia first diverges from the equilib- 
rium curve at the N2-NH3 interconversion quench point, 
which is at relatively high temperatures and pressures 
due to the strong N=N triple bond. However, there is 
a pressure and temperature dependence to the different 
chemical pathways that allow NH3 -o- N2 interconver- 
sion, and different reactions begin to dominate at differ- 
ent pressure levels. More importantly, the temperatures 
in the HD 189733b models do not drop off significantly 
between the "initial" NH3 quench point near '^200 bar 



and the "final" quench point at the few-bar level, allow- 
ing the interconversion reactions to continue to operate 
at a reduced but still effective rate in the intervening re- 
gion. The final NH3 quench occurs when the thermal 
gradients become large again and the temperatures be- 
come low enough that the NH3-N2 interconversion reac- 
tions cease to become important. Ammonia quenching 
then operates over an extended range on HD 189733b, 
as marked by the two quench-point arrows in the HD 
189733b NH3 panel in Fig. [31 The thermal gradients are 
larger in the region where NH3 initially quenches on HD 
2G9458b, allowing chemical conversion reactions to shut 
down more abruptly, with a better-defined single quench 
level for NH3 in the HD 209458b models. 

The details of the particular quenching mechanisms, 
and their sensitivity to temperature and transport as- 
sumptions, will be discussed below. At this point, we 
simply want to illustrate that transport-induced quench- 
ing can be an important process on extrasolar giant plan- 
ets that can allow disequilibrium species to be carried 
to relatively high altitudes to potentially affect spec- 
tral behavior. Note that transport-induced quenching 
is much more important on the cooler HD 189733b than 
on HD 209458b. Ahhough Fig. [3] shows that quenching 
does occur in our dayside HD 209458b kinetics/transport 
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models, the effectiveness of transport-induced quenching 
on HD 209458b is limited by the presence of a strong 
stratospheric thermal inversion — as temperatures be- 
come hot enough in the stratosphere, thermochemical 
kinetics drives the composition back toward equilibrium. 
When the inversion is weak or absent on HD 209458b 
[i.e., at the terminators or nightside), transport-induced 
quenching becomes more important. 

The stratospheric thermal inversion also affects the 
height to which the molecules can be carried. The models 
for both planets in Fig. [3] have the same eddy diffusion 
coefficient, yet a comparison of the dotted-line models 
for both planets shows that the species reach different 
pressure levels on each planet. The difference is only 
partly due to the higher temperatures creating a lower 
homopause level on HD 209458b (where the term "ho- 
mopause" refers to the pressure level at which the molec- 
ular and eddy diffusion coefficient are equal; above the 
homopause, molecular diffusion acts to limit the abun- 
dance of heavier species, whereas below the homopause 
the atmosphere can be well mixed). This point is demon- 
strated by a close comparison of CO and N2 on each 
planet, as these thermochemically stable species reach 
their homopause levels, yet they are carried to only 
slightly different pressures on each planet. The main 
cause of the limited penetration of molecular species to 
high altitudes on dayside HD 209458b is the high strato- 
spheric temperatures — species like CH4 and NH3 can- 
not survive the thermochemical kinetics at these temper- 
atures and are converted in to atomic C, N, and H, as 
well as to the more stable molecular species CO and N2. 

Another interesting point to note from Fig. [3] is that 
some molecules like CO, H2O, N2, and CO2 are rela- 
tively unaffected by transport-induced quenching on ei- 
ther planet. Although CO, H2O, and N2 technically 
also quench in the thermochemical kinetics and trans- 
port models, they are already the dominant species at 
their quench points, and the quenching does not greatly 
affect their abundances. As we will show below, this 
result is somewhat dependent on assumptions about 
the eddy diffusion coefficient and temperature structure. 
On HD 189733b, CO2 quenches eventually (see where 
dotted curve diverges from dashed curve), but only at 
very high altitudes, and its abundance follows equilib- 
rium at the pressure levels at which th e observations 
are th e most sensitive (^I-IO""* bar, see iFortnev et all 
l2008al) . Other species like hydrocarbons, nitriles, and 
most nitrogen-bearing species are strongly affected by 
transport-induced quenching. 

A comparison of the solid curves with the dotted curves 
in Fig. [3] shows the influence of photochemistry. When 
ultraviolet photons are present, molecules like H2O, CO, 
N2, and NH3 are photolyzed, initiating a slew of photo- 
chemical reactions th at alter the chemis try and composi- 
tion. As discussed bv lLiang et al.l ()2003[ ). one main result 
of the photochemistry in the atmospheres of close-in ex- 
trasolar giant planets is the production of huge quantities 
of atomic hydrogen from photolysis of H2O, followed by 
catalytic destruction of H2 from the reaction of molecular 
hydrogen with the photolytic product OH. Other similar 
catalytic cycles are initiated deeper in the atmosphere 
from N H3 photolysis and fr om reaction of H with CH4 
(see alsol Zahnle et al.ll2011[ ). 

Because many molecules are vulnerable to attack by 



atomic hydrogen at these temperatures, the production 
of atomic H sets the stage for the subsequent photo- 
chemistry. Some molecules like H2O, N2, and CO are 
efficiently recycled and are relatively stable against pho- 
tochemistry. As can be seen from Fig. [31 photochemistry 
only affects the profiles of these species at very high alti- 
tudes, where visible and infrared transit and eclipse ob- 
servations have little sensitivity. Other molecules like 
CH4 and NH3 have weaker bonds and/or are more re- 
active and less efficiently recycled. These molecules can 
be removed from portions of the visible atmosphere, to 
be replaced by photochemical products. The big win- 
ners in terms of photochemistry are atomic species such 
as H (by far the dominant photochemical product), C, 
O, and N; some radicals like CH3, NH2, and OH; un- 
saturated hydrocarbons like C2H2; and some nitriles like 
HCN. Note that photochemistry (and thermochemistry, 
if the atmosphere is warm) can keep molecular species 
from reaching their homopause levels, as species can be- 
come kinetically destroyed faster than eddy diffusion can 
transport them to high altitudes. Below we discuss the 
kinetic mechanisms in more detail. 

3.2. Quench kinetics 

As shown in Fig. |3] and discussed above, transport- 
induced quenching is an important process in the atmo- 
spheres of extrasolar giant planets. For relatively warm 
planets like HD 189733b and HD 209458b, the thermal 
structure lies within the CO and N2 stability fields at 
the temperatures and pressures where quenching occurs, 
so that the less abundant CH4 and NH3 are strongly af- 
fected by quenching. This situation differs from that of 
cooler planets like GJ 436b (which is in the CH4 and N2 
stability fields, so that quenching of CO and NH3 can 
be important) or from cold planets like Jupiter (which is 
in the CH4 and NH3 stability fields, so that quenching 
of CO and N2 are important) or from brown dwarfs like 
Gliese 229b whose thermal structure crosses the CO = 
CH4 and N2 — NH3 boundaries in such a location that 
quenching of all four molecules — CO, CH4, N2, and 
NH3 — could be important. 

The kinetic pathways leading to the quenching of 
CO-CH4 interconversion, N2-NH3 interconversion, and 
NH3-HCN-CH 4 interc onve rsion have been di scussed in 
iVisscher et all (|2010bf ) and lMoses et"all (|2010f ). We have 
made several updates to the rate coefficients used in those 
studies, including adopting the CH3 -|- OH reaction path- 
ways and rate coefficients of IJasper et al.l (j2007D , and 
updating the rate coefficients for th e nitrogen reactions 
based on jClippcnstcin ct all (|2009f) and Klippenstem] 
(|2010D . Because of uncertainties in the kinetic rate coef- 
ficients for certain key r eactions in the CO O CH4 inter- 
conversion scheme fsee IVisscher et al]|2010bl ). we have 
taken a more detailed look at the CH3OH -I- H reaction 
and its various pathways for this paper. In particular, 
we have used ab initio transition-state theory to predict 
the rate coefficients for the three primary channels of the 
reaction of H with CH3OH: 

R864 CH3OH + H^ CH2OH + H2 
R862 CH3OH + H^ CH3O + H2 
R860 CH3OH + H^ CH3 + H2O. 
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We used the Q CISD(T)/cc-py TZ method (e.g. 

iPople et al.l Il987t iDunning etall 119890 to predict 
the rovibrational properties of the stationary points. 
Higher-level energies were obtained from basis-set 
extrapolation of exphcit QCISD(T) calculations 
for the cc-pVTZ and cc-pVQZ basis sets. These 
QCISD(T)/CBS//QCISD(T)/cc-pVTZ calculations 
were performed with the spin-restricted formalism using 
the MOLPRO software packageQ Zero-point corrected 
barrier heights of 11.0, 13.4, and 24.2 kcal/mol were 
obtained for channels R864, R862, and R860, respec- 
tively. Most of the modes were treated as rigid-rotor 
harmonic oscillators, but hindered-rotor torsional and 
asymmetric Eckart-tunneling corrections were included 
as appropriate. 

The predicted rate coefficients for reactions R864- 
R860 are well reproduced over the 500 to 2500 
K region by the modified Arrhenius expressions, 
ks64 = 1.09 X 10-19 T2-^28 gxp(_2240/T), ks62 
= 6.82 X 10-20 2^2.658 gxp(_4g43/2.)^ ^^^^ ^ 

4.91 X 10-19 exp(-10380/T), cm^ molecule-^ 

s"!. These predictions are expected to be accurate 
to within a factor of two. For channels R864 and 
R862, the predictions are wit hin a factor of 1 .35 o f 
the G2-based estimates from iJodkowski et al.l ()1999l ) 
but have lower uncertainties due to the use of higher 
levels of elect ronic-structure t heory . The more recent 
predictions of iCarvalho et al.l ()2008[ ) are closely related 
to the present ones, being only about 0.8 times lower. 
The only prior study of channel R860 app ears to be 
the G2-lDased study of iLendvav et al.l (|1997[ ). who did 
not provide explicit expressions for this rate coefficient. 
The absence of both theoretical and experimental data 
for channel R860 was the motivating factor for the 
present ab in i tio tra nsition-state theory analysis, as 
iVisscher et al.l (|2010bl ) have suggested that this pathway 
could be important in the overall interconversion scheme 
between CH4 and CO on giant planets. 

The above ab initio calculations indicate, however, 
that the reaction R860 (CH3OH + 11^ CH3 + H2O) is 
too slow to participate in the dominant CH4-CO quench- 
ing scheme in our HD 189733b and HD 209458b models. 
Instead, we find that the dominant reaction scheme con- 
trolling CH4 — > CO conversion and methane quenching 
depends on the eddy diffusion coefficient and pressure- 
temperature conditions at the quench point; the domi- 
nant scheme can even bypass methanol (CH3OH) alto- 
gether in some cases. For smaller eddy diffusion coeffi- 
cients, the quench level is at lower pressures (higher al- 
titudes), and the dominant CH4-CO conversion scheme 
becomes 



H 
H 



CH4 
H2O 



CH3 + H2 
OH + H2 



^ MOLPRO is a package of ab initio programs written by H.- 
J. Werner, P. J. Knowles, F. R. Manby, M. Schiitz, P. Celani, 
G. Knizia, T. Korona, R. Lindh, A. Mitrushenkov, G. Rauhut, T. 
B. Adler, R. D. Amos, A. Bernhardsson, A. Berning, D. L. Cooper, 
M. J. O. Deegan, A. J. Dobbyn, F. Eckert, E. GoU, C. Hampel, 
A. Hesselmann, G. Hetzer, T. Hrenar, G. Jansen, C. Koppl, Y. Liu, 
A. W. Lloyd, R. A. Mata, A. J. May, S. J. McNicholas, W. Meyer, 
M. E. Mura, A. NicklaB, P. Palmieri, K. Pfiiiger, R. Pitzer, M. Rei- 
her, T. Shiozaki, H. StoU, A. J. Stone, R. Tarroni, T. Thorsteins- 
son, M. Wang, and A. Wolf. 



OH + CH3 
CH2OH + M 
H2CO + H 
HCO + M 



CH2OH + H 
H + H2CO + 
HCO + H2 
H + CO + M 



M 



Net : CH4 + H2O ^ CO + 3H2, 



(2) 



where M represents any third body. For larger eddy dif- 
fusion coefficients, the quench level is at higher pressures 
(lower altitudes), and the dominant CH4-CO conversion 
scheme can become 



H 
H 

OH + CH3 
CH3OH 
CH2OH 
H2CO 
HCO 
H2 



CH4 — !• CH3 + H2 

H2O OH H2 
M-^ CH3OH + M 
H ^ CH2OH + H2 
M ^ H + H2CO + M 
H ^ HCO + H2 
M ^ H + CO + M 
M ^ 2 H + M 



Net : CH4 + H2O ^ CO 



3H, 



(3) 



The rate-limiting step for the CH4 CO conversion 
process is the slowest reaction in the fastest overall con- 
version scheme. For the above scheme (2), the rate- 
limiting reaction is OH + CH3 CH2OH -I- H (reac- 
tion number R670 in our list), and for the above scheme 
(3), the rate-limiting reaction is OH -I- CH3 -I- M — >■ 
CH3OH + M (reaction R674). The relative importance 
of schemes (2) and (3) depends on which of these OH -f 
CH3 pathways is faster for the pressure-temperature con- 
ditions at the quench level for the planet under consid- 
eration. We adopt rate-coefficient expressi ons for these 
two p otential rate-limiting reactions from iJasper et al.l 
(|2007f) . who provide a thorough description of the pres- 
sure and temperature dependence. Both schemes typi- 
cally contribute some nontrivial fraction of the rate of 
interconversion in our models, so considering both path- 
ways in a combined rate- limiting step can furnish a good 
description of the methane mo le-fraction quenching be - 
havior on transiting exoplanets (jVisscher fc Mosej2011|) . 

Note that although the C-H-0 system has been well 
studied due to applications in combusti on chemistry 
and terrestrial atmospheric chemistry (e.g.. lBaulch et aP 
l2005t I Atkinson et 311120061) . the kinetics of CO O CH4 
interconversion in reducing environments remains some- 
what uncertain. The rate coefficients for the above rate- 
limiting reactions are uncertain by perhaps a factor of 3, 
and other alternative schemes may be competing if other 
rate coefficients in our full list are in error. The methane 
mole fraction on HD 189733b and HD 209458b will de- 
pend critically on the kinetic rate coefficients adopted in 
the model, particularly for the reactions that form C-0 
bonds from species with 0-H and C-H bonds. The abun- 
dance of quenched species is always highly dependent on 
the assumed kinetic rate coefficients, as well as on the 
assumed strength of atmospheric mixing. 

Despite these uncertainties, the OH + CH3 + M 
^ CH3OH + M and the OH + CH3 ^ CH2OH 
+ H reactions are much more likely to be the rate- 
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limiting step in the conversion of CH4 to CO on exo- 
planets than the reaction OH + CH3 ^-112 + H2CO 
(i.e., the reverse of the CO quench re action origi- 
naUy proposed by iPrinn fc Barshavl Il977| ) , as the lat- 
ter reaction is a few orders of magnitude too slow to 
be a likely player in any dominant con version s cheme 
see iDean & Westmoreland 1987: Yu ng et al.l 119881: 



nii^et 

Griffith & Yell e ,1999^ .Bezard et al., .2003 iJasper et al.l 
20071: IVisschcr et al.l l2010b[) . Note also that the dom 



inance of CH2OH over CH3O as a product in the re- 
action of H with CH3OH and in the reaction of OH 
with CH3, combined with the rapid thermal decom- 
position of CH2OH, suggests that the dominant path- 
way will involve CH2OH as an intermediate, not CH3O, 
which was suggested or adopted in p revious CO or CH4 
quen c hing studies fe.g.. lYung et alT l988: Bezard ctaD 
2OOI; 'CooTJer & Showman' '2006'; 'Visscher et aL| 12010^ 
Line et al..,2010: Madhusudhan fc Seager.,201lll 



The transpo rt time scale arguments of 

IPrinn fc Barshavl ()1977[) can be used to predict the 
quenched abundance o f CH4 on exoplanets (see also the 
updates proposed by IVisscher fc Moses! 120 111) . With 
this procedure, the quench point is located where the 
chemical kinetic conversion time scale Tchem is equal 
to the vertical transport time scale Tdyn- 

Then, the 

mole fraction of the quenched species above the quench 
level is simply the equilibrium mole fraction of the 
constituent at that quench point. In our models, the 
chemical kinetic conversion time scale for CH4 is 



[CH4 



d[CH4]/di 



[CH4 



fc670[CH3] [OH] + fc674[CH3] [OH] [M] 



(4) 



where [X] represents the equilibrium concentration of 
species X in cm"'^, kgro is the rate coefficient for reaction 
R670 in cm^ s~^, and kfi74 is the rate c oefficient for R670 
in cm s~ ^ See IJasper et al.l (|2007D . IVisscher fc MosesI 
(|2011|) . and Supplementary Table SI for more details on 
kgyo and k674. 

The transport time scale is Tdyn — L"^ /K^z, where L 
is an effective length scale over which dy namical mixing 
operates. As discussed bv iSmithI ()1998l) . L is generally 
some fraction of the scale height H . For quenching of 
CO on cold planets like Jupiter, L w 0.1-ff, whereas for 
quenching of CH4 o n close-in transiting hot Jupiters, L 
is closer to Q.hH (see IVisscher fc Mosesll201ll) . When we 
use the procedure suggested bv lSmithI (|1998D to derive L 
and when we calculate Tchem as described in Eq. (4) for 
our rate-limiting step, we derive a quenched mole frac- 
tion within ^^15% of our kinetics/transport model results 
for HD 1897 33b and HD 2 9458b . Therefore, as is also 
discus s ed bv iBezard et all (I2002D . iCooper fc ShowmanI 
(l2006l). IVisscher et all (l2010bD . iMoses et al.l (|2010[). and 
IVisscher fc MosesI ( 2011| ) the tra nsport time-s cale ap- 
proach is valid provided that the iSmithI (|1998l ) formal- 
ism for Tdyn is used and provided that the assumed ki- 
netic rate-limiting step is a reasonable one. We cau- 
tion that the time-scale approach has been misused fre- 
quently in exoplanet and brown-dwarf literature, where 
errors in the calculation of the rate coefficients for 
the rate-limiting steps have arisen due to incorrect re- 



versals of three -body reactions (jCriffith fc Yelld 119991 : 
iLine et aI]|2010D . or where numerous other misconcep- 
tions about rate-limiting reactions and their rate coef- 
ficients and transport time scales are ub i quitous (see 
Bezard et al. 2002, Visscher et aLj I2010bl . IMoses et all 



201a or IVisscher fc Moses 201irfor a more detailed 



discussion of the problem). Some of these problems 
have serious implications with respect to conclusions 
about the strengt h of atmospheric mixing, particularly 
on br own dwarfs (jTrover et al.ll2007t IVisscher fc MosesI 
120 ill ), where K^z is likely several orders of magni- 
tude greater than the sluggish Kzz ~ 10^-10^ cm^ 
s~^ that has often been re p orted ( e.g.. IGriffith fc Yelld 
19991 ISaumon et al.l I200I I200I ILeggett et al.l 120071: 



Geballe et al. I 120090 . but also in relation to the im- 



plied strength of atmospheric mixing on exoplanets (e.g. 
iMadhusudhan fc Seageill201lD . 

The time-scale approach fails for species that have 
more than one quench point, for species that have an 
extended quench region due to Tchem remaining close to 
Tdyn over an extended pressure region, or for species that 
otherwise exhibit complicated kinetic behavior due to 
quenching of other key species. As an example, NH3 
in our HD 189733b models quenches in a region of the 
atmosphere in which the temperature does not fall off 
sharply with height. The NH3 mole fraction departs from 
the equilibrium curve when Tdyn drops below Tchem and 
conversion between NH3 and N2 first becomes inhibited; 
however, due to the nearly isothermal behavior of the 
atmosphere above the quench region, Tchem remains rel- 
atively close to Tdyn over an extended pressure range, and 
NH3 does not fully quench until the temperatures begin 
to fall off sharply again. Therefore, the time-constant ar- 
guments provide a good prescription for where the NH3 
leaves equilibrium but not for the final quenched mole 
fraction. 

Several mechanisms closely compete as the dominant 
mechanism involved in NH3 — >■ N2 conversion, and again, 
the dominance depends on pressure-temperature con- 
ditions near the quench point. In our HD 209458b 
models, where the atmosphere is relatively warm and 
NH3 quenches at relatively low pressures, the dominant 
scheme is 



2 (NH3 + H 
NH2 + H 

NH + NH2 
N2H2 + H 
NNH + M 



NH2 + H2) 
NH + H2 
N2H2 + H 
NNH + H2 

N2 + H + M 



H2-t-M^2H-fM 
Net : 2NH3 N2 -I- 3H2. 



(5) 



On the cooler HD 189733b, NH3 first quenches down 
deep at high pressures and slightly lower temperatures 
than at the quench point on HD 209458b, through the 
dominant scheme 

NH3 + H ^ NH2 + H2 
NH2 -f NH3 N2H3 + H2 
N2H3 -f M ^ N2H2 + H + M 
N2H2 + H ^ NNH + H2 
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NNH + M^N2+H + M 
Net : 2NH3 ^ N2 + 3H2. 

(6) 

However, this scheme (6) becomes less important with 
decreasing pressure, and scheme (5) is responsible for the 
final quench point on HD 189733b. The rate-limiting 
step in schen ie (5) is NH + NHp -» N2H2 + H, and 
we adopt the iKlippenstein et ahl (|2009D theoretical rate 
coefficient expression of 7.07 x 10"!° T-°-272 exp(39/T) 
for this reaction. In scheme (6), the rate-limiting step 
is N2H3 -I- M ^ N2H2 + H + M, where we take 
th e rate coefficient from the t heoret ical calculations 
of iDean fc Westmorelandl ()1987D and iHwang fc Mebell 
(12003!) (see Supplementary Table SI). Note that these 
two reaction schemes — and all other potential NH3-N2 
conversion schemes — are somewhat speculative and rely 
on rate coefficients derived from theoretical calculations 
rather than laboratory measurements. As a whole, the 
kinetics of nitrogen compounds is less well understood 
than that of the C-H-0 system (sec Moses et al. 2010 for 
further discussion). 

After conversion of NH3 to N2 is first quenched (at 
-1600 K in the HD 189733b models, and at >1700 K in 
the HD 209458b models), the disequilibrium abundance 
of NH3 affects the abundances of other nitrogen-bearing 
constituents, often at altitudes well below their ultimate 
quench points. HCN is a prime example. The HCN mole 
fraction departs from equilibrium at the first NH3-N2 
quench point, but HCN remains in a pseudo-equilibrium 
with NH3, CH4, and H2 until its final quench point at 
higher altitudes. The dominant scheme involved in NH3 
— > HCN conversion in our cooler HD 189733b models is 

NH3 + H ^ NH2 + H2 
CH4 + H ^ CH3 + H2 
NH2 + CH3 + M-¥ CH3NH2 + M 
CH3NH2 -h H ^ CH2NH2 + H2 

CH2NH2 -> CH2NH -I- H 
CH2NH + H ^ H2CN -f H2 
H2CN + M ^ HCN + H + M 

H2+M-^2H-HM 

Net : NH3 + CH4 HCN + 3H2. 

(7) 

This scheme is the reverse of the dominant HCN 
NH3 conv ersion schenie der ived for Jupiter's deep tro- 
posphere (|Moses et al.ll20ldl ). although the rate-limiting 
step in the above scheme (7) for our exoplanet model- 
ing is CH3NH2 4- H - > CH2NH2 + H2. We adopt the 
IDean fc Bozzellil (|2000[ ) theoretically-derived rate coef- 
ficient of 9.3 X 10-i'^Ti-5exp(-2750/r) cm^ s^^ for 
this rate-limiting reaction. On the warmer HD 209458b, 
HCN exchanges carbon with CO, as well as with CH4, 
and HCN is quenched when the following scheme be- 
comes ineffective: 

NH3 -I- H NH2 + H2 
CO + NH2 HNCO + H 
HNCO H HCN + OH 
OH H2 ^ H2O + H 



Net : NH3 -f CO ^ HCN + H2O. 

(8) 

Some molecules like C2H2 quench at very high alti- 
tudes or are never truly quenched, as kinetic processes 
continue to affect their abundances throughout the at- 
mospheric column, even when photolyzing radiation is 
absent. The dominant mechanism controlling the inter- 
change of C2H2 and CH4 in our models is the following 
scheme or its reverse: 

C2H2 -f H + M^ C2H3 + M 
C2H3 + H2 — > C2H4 + H 
C2H4 + H + M^ C2H5 + M 
C2H5 + H2 ^ C2H6 + H 
C2H6 + M ^ 2 CH3 + M 
2(CH3 + H2 ^ CH4 + H) 

2 H + M -» H2 -f M 

Net : C2H2 -I- 3 H2 ^ 2 CH4 . 

(9) 

The rate-limiting step in this scheme is either C2H6 ther- 
mal decomposition or H2 -t- C2H5 -> C2H6 -I- H, depend- 
ing on local pressure-temperature conditions. 

In all, transport-induced quenching can have a signif- 
icant effect on the predicted abundances of species like 
CH4, NH3, HCN, and C2H2 on hot Jupiters. Given the 
overall uncertainties in the thermal structure, eddy dif- 
fusion coefficients, kinetic rate coefficients, and effects 
from horizontal transport, the quantitative predictions 
from our kinetics and transport models should not be 
taken too seriously. However, the general behavior in 
the models should be heeded — the column abundance 
of these species will likely be enhanced over equilibrium 
predictions (and in some cases significantly enhanced), 
and transport-induced quenching is more significant for 
cooler planets or atmospheric regions. 

We now move on to discuss the effects of the photo- 
chemistry that is initiated when ultraviolet photons are 
absorbed by atmospheric constituents. 

3.3. Photochemistry of oxygen species 

Figure |4] shows vertical profiles of the major oxygen- 
bearing species in our nominal thermo/photochemical ki- 
netics and transport models. Photolysis of the major 
parent molecules CO and H2O drives much of the oxy- 
gen chemistry. Water photolysis in particular is a key 
component of exoplanet photochemistry, as it sets up a 
catalytic cycle that destroys H2 and produces H: 

H2O + /li/ ^ H + OH 
OH + H2 ^ H2O + H 
Net : H2 ^ 2 H, 

(10) 

where hiy represents an ultraviolet photon. Water is re- 
cycled with this efficient scheme, but H2 is destroyed, 
resulting in a huge net production rate for atomic H, 
which then replaces H2 as the dominant atmospheric 
constituent at high altitudes (see Fig. 0]). This transi- 
tion from H2 to H likely signals a transition into the 
high-temperature thermosphere on hot Jupiters, as no 
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potential molecular coolants survive when atomic H be- 
com es the dominant con stituent at high altitudes (see 
also lGarcia Mufio^l2007D . Atomic species therefore are 
expected to dominate in the t hermospher e s of these close- 
in transiting planets (e.g. , lYelld 120041 iGarcfa Munod 
I2007t iKoskinen et all l2010l ). which might help explain 
the lack of obvious H2 emission in far-ultraviolet day - 
glow observations of HD 209458b (jFrance et al.ll201oh . 
From a column-integrated sense, scheme (10) is signif- 
icantly more important as a source of atomic H than 
direct H2 photodissociation. Atomic H flows downward 
from its production region and remains at abundances 
greatly in excess of equilibrium until in reaches altitudes 
marked by its mole-fraction minimum shown in Fig. |4l 
near 0.01-0.1 bar on both planets, at which point atomic 
H again follows its equilibrium curve to lower altitudes. 

Photolysis of CO to produce C + O and photolysis of 
H2O to produce 2 H + O or 0(^D) + H2 also contribute 
to H2 destruction through schemes such as 

CO + hv^C + 
+ H2 -> OH + H 

OH + H2 ^ H2O +H 

Net : CO + 2H2 ^ C + H2O + 2H 

and 

H2O + /li^ ^ 2H + 
+ H2 ^ OH + H 
OH + H2 ^ H2O + H 
Net : 2H2 ^ 4H 

and 

H2O + hv ^Yl2 + 0(^D) 
0(iD) -h H2 ^ OH -l-H 
OH -H H2 ^ H2O + H 
Net : H2 ^ 2 H. 

Reactions of atomic H with H2O and OH can partially 
recycle the H2 , but key reactions like 2 H -|- M ^> H2 
do not operate effectively at low pressures, and the net 
result is that the loss rate for H2 exceeds the production 
rate at high altitudes. The H atoms produced in these 
mechanisms diffuse both upward and downward, where 
they strongly affect the abundance of other species. 

Although H2O photolysis operates effectively down to 
millibar regions in these models and the atomic H pro- 
duced from the photolysis continues to attack the H2O, 
water is not permanently removed from this atmospheric 
region because of the efficient recycling by the reaction 
OH -I- H2 — >■ H2O -I- H. Note that rapid replenishment 
by diffusion from other altitudes can also play an impor- 
tant role in maintaining H2O for our models with cooler 
stratospheric temperature profiles (see also iLine et aD 
|2010( ). As an example, water is in approximate chem- 
ical equilibrium at pressures greater than '^O.l bar in 
the HD 189733b dayside-average model shown in Fig. |31 
however, the transport time scales are shorter than the 
net photochemical loss time scales for H2O in the region 



from '^0.1 bar to a few /ibar, so that atmospheric mixing 
"quenches" the H2O mole fraction and keeps it roughly 
constant through this region. Since the H2O mole frac- 
tion quenches at a value similar to the equilibrium value 
and since equilibrium would keep the H2O mole frac- 
tion constant in this region anyway, the water profile for 
the thermo/photochemical kinetics and transport model 
is quite similar to the thermochemical equilibrium pro- 
file for our HD 189733b models. For our much warmer 
dayside HD 209458b model, equilibrium chemistry dom- 
inates the H2O profile over a much more expanded alti- 
tude range — transport is only important at high alti- 
tudes, and only if the eddy diffusion coefficient is large. 
This effect can be seen from detailed comparisons of the 
H2O profile from the HD 209458b dayside-average model 
in Fig.[3]in which it is assumed that Kzz = 10^ cm^ s~^, 
in contrast to the H2O profile in Fig. |4] from an otherwise 
similar HD 209458b dayside-average model in which it is 
assumed that mixing is more vigorous at high altitudes 
(see the nominal K^z profile from Fig. [1]). The Kzz — 
10^ cm^ s~^ model has the H2O profile closely following 
equilibrium, whereas our higher Kzz nominal model has 
water being carried to higher altitudes before its even- 
tual chemical destruction. However, this additional H2O 
from transport processes does not add much to the col- 
umn abundance of H2O on the planet. The net result 
on both the cooler and warmer classes of planets is an 
H2O profile that closely follows the equilibrium profile 
through the observable regions of the atmosphere. 

Carbon monoxide also survives in the atmospheres of 
HD 189733b and HD 209458b despite rapid loss from H 
-h CO M ^ HCO + M, from OH -I- CO ^ CO2 -f 
H, and from photolysis by EUV photons that continues 
down to mbar levels. The first two reactions dominate 
the column-integrated CO loss in the stratosphere; the 
HCO produced in the first reaction efficiently reacts with 
H to reform CO, and the CO2 produced in the second 
reaction efficiently reacts with H to reform CO, so carbon 
monoxide is readily recycled on both HD 189733b and 
HD 209458b. Even in the case of CO photolysis, some 
recycling pathways such as the following can operate: 

CO + hv^G + & 
C + H2 -> CH -f H 
CH -f H2O ^ H2CO -t-H 
H2CO + H ^ HCO + H2 
HCO -H H ^ CO + H2 
Net : H2O ^ + H2, 

allowing CO to persist. However, these recycling schemes 
are not 100% effective, and some of the CO photodestruc- 
tion leads to the production of other species. Atomic 
carbon and oxygen released from CO photolysis at high 
altitudes can remain as C and O, for example, and some 
of the O can lead to H2O production through reactions 
like O H2 ^ OH -h H followed by OH + H2 ^> H2O 
-|- H, while some of the C can go tow ard producing hy- 
drocarbons and HCN (see below, a nd lLiang et al.ll2004l : 
ILine et al.|[20lol: iZahnle et all 1201 1[ ). As with H2O, the 
CO profiles are expected to closely follow the equilibrium 
profiles, except at high altitudes, where the CO will be 
converted to C, O, and H2O (i.e., for cooler-temperature 
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Fig. 4. — Mole-fraction profiles for several oxygen-bearing species in our models for HD 189733b (left) and HD 209458b (right) for 
dayside-average (top) and terminator-average (bottom) conditions, assuming the nominal K^z profiles shown in Fig. [T] A color version of 
this figure is available in the online journal. 



models) or where high eddy diffusion coefficients may 
allow it to be transported to higher altitudes than equi- 
librium models predict (i.e., for our nominal dayside HD 
209458b model). 

The photochemistry of CO2 is not particularly inter- 
esting, as CO2 remains in equilibrium with CO and H2O 
throughout most of the atmosphere, via the reaction 
scheme: 

H2O + H OH +H2 
OH + CO o CO2 + H 



Net : H2O + CO o CO2 



H2 



Some additional CO2 in excess of equilibrium is pro- 
duced at high altitudes from CO photochemistry on HD 
189733b and/or from transport on HD 209458b, but the 
CO2 profiles on both planets remain close to equilibrium 
predictions. 

Aside from atomic H, few photochemical products 
build up from the chemistry of oxygen-bearing species 
in our hot- Jupiter models. Some O, OH, NO, and even 
O2 and HNCO molecules are produced at high altitudes 
and diffuse downward; however, the photochemical life- 
times of these species are quite short, and they are not 



likely to achieve column abundances great enough to be 
detectable at visible or infrared wavelengths with current 
technologies unless the planet's metallicity is enhanced 
greatly over solar (see Fig. |4] and Table [Ij . Atomic O 
is produced from the reaction of H with OH to form 
O -I- H2, with a lesser contributions from H2O, OH, 
and CO photolysis. Atomic O is lost predominantly 
from O + H2 — > OH -I- H. Because atomic O is not 
readily destroyed by H, it can be transported up to 
its homopause level in our HD 209458b models, unlike 
the situation for molecular species. Some fraction of 
the O column will then make it into the thermosphere 
on planets with vigorous vertical mixing, where high 
thermospheric temperatures and hydrodynamic winds 
or other processes may then distribute the ato mic O 
over a large radial distance ( Lccavclier des Etangs et alj 
2004 lYelle|l2004t iTian et al.1,2005; .G arcia Mu iioa 12007 



Penz et al.ll2008tlMurrav-Clav et all2009: La mmer et al.l 



20091 ) ■ which will increase the likelihood of its detection 



at ul t raviolet wavelengths (e.g., see [V idal-Madiar et al] 
[2OOI lBen-Jaffel|[20ToriKoskinen et all I2010f ). "'Even in 
cases where K^^ is small in the upper atmosphere, the 
dynamical wind from the hydrodynamic escape — which 
is not included in our models — may then dominate over 
molecular diffusion in controlling the behavior of heavy 
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atmospheric constituents, so that the concept of a ho- 
mopause may not actually be appropriate for close-in 
transiting planets. The dominant OH production mecha- 
nism is H + H2O — > OH -I- H2, with the atomic H coming 
from the H2 catalytic destruction scheme (10) described 
above. The reverse of this reaction is its dominant loss 
mechanism. In general, the cooler the stratosphere, the 
larger the increase in the OH abundance over equilibrium 
predictions. 

The main production mechanism for NO is N -I- OH 
NO -I- H. The OH derives from water, and the N can 
come from either N2 photolysis or from NH3 via schemes 
such as the following: 

H2O + hiy ^ OU + U 
NH3 + H ^ NH2 + H2 
NH2 H ^ NH + H2 
NH -I- H ^ N + H2 

N OH NO H 

Net : H2O + NH3 H -J- NO 3H2. 



The double-peak appearance in the NO mole fraction is 
due to the dominance of the N2 photolysis pathway in the 
upper stratosphere and the NH3 pathway in the middle 
stratosphere. The NO is lost through reactions with N to 
form N2 + O, with C to form N CO and/or CN + O, 
through photolysis to form N -|- O, and through reaction 
with atomic H. Molecular oxygen is produced through 
the reaction O + OH — > O2 -I- H, and is lost through the 
reverse of this reaction. 

As a general rule of thumb, photochemistry in 
these warm, highly irradiated stratospheres favors small 
molecules and/or molecules with strong bonds. Our re- 
sults are co nsistent with othe r exop l anet photochem- 
ical model s fLiang.^LalJ I2OOI [2OOI iLine et al.l [20T0I : 
iZahnle et~a l. 2011) in that we see little photochemical 
production of molecules like H2CO, CH3OH, H2CCO, or 
CH3CHO. 

3.4. Photochemistry of carbon species 

The photochemistry of carbon species is more inter- 
esting than that of oxygen species, as the primary par- 
ent molecule CH4 has weaker bonds. Several carbon- 
bearing photochemical products can build up in the at- 
mospheres of extrasolar giant planets. Figure [5] shows 
vertical profiles of the major carbon species in our nom- 
inal thermo/photochemical kinetics and transport mod- 
els. Carbon photochemistry is initiated by CO photol- 
ysis at high altitudes and by hydrogen abstraction from 
methane at lower altitudes. At pressures less than ~1 
/ibar, much of the C that is released from CO photolysis 
remains as atomic C, although some acetylene and other 
hydrocarbon production can occur through schemes such 
as 

2 {CO + hiy ^ C + 0) 
2(C + H2 ^ CH +H) 
CH + H2 ^ ''CH2 + H 
CH + 3CH2 ^ C2H2 + H 
2 ( 2 H + M ^ H2 M ) 



Net : 2C0 + H2 ^ C2H2 + 2 0. 

(11) 

At greater pressures (and in particular where CH4 or 
CH3 become more abundant than atomic C), acetylene 
and other hydrocarbons are produced from methane de- 
struction through pathways such as 

2(H20 + hi^ ^ U + OE ) 
2 (OH + H2 ^ H2O +H) 
2 ( CH4 + H ^ CH3 + H2 ) 
CH3 + H — )■ ^CH2 + H2 
"'^CH2 + H2 — > ■^CH2 + H2 
3CH2 + CH3 ^ C2H4 + H 
C2H4 + H — > C2H3 + H2 
C2H3 + H — > C2H2 + H2 
Net : 2CH4 ^ C2H2 + 3H2. 

(12) 

Note the importance of H2O photolysis in providing the 
atomic H used to destroy the CH4. These pathways are 
efficient enough that C2H2 becomes the dominant hydro- 
carbon at high altitudes under certain conditions in our 
models (see Fig. [5]). 

Acetylene is more effectively hydrogenated with in- 
creasing pressure, such that the following scheme acts 
to remove the C2H2 and recycle the methane at middle- 
stratospheric pressures: 

C2H2 + H + M ^ C2H3 + M 
C2H3 + H2 — > C2H4 + H 
C2H4 + H + M ^ C2H5 + M 
C2H5 + H ^ 2 CH3 
2 ( CH3 + H2 ^ CH4 + H ) 
Net : C2H2 + 3E2 ^ 2CH4. 

(13) 

This scheme, rather than photolysis, is responsible for 
removing C2H2 from the middle and lower stratosphere, 
as the C2H2 photolysis products tend to simply recycle 

the C2H2. 

Other C2H2; hydrocarbons are less able to survive the 
large background H abundance. Ethane is produced in 
the stratosphere through 2 CH3 + M ^ C2H6 -I- M and 
through C2H5 -I- H2 — > C2H6 + H. Ethane is lost through 
hydrogen abstraction by atomic H (C2H6 -I- H ^ C2H5 
+ H2), followed hy a + C2H5 2CH3, and eventual 
methane production. The column abundance of photo- 
chemically produced ethane is quite small in our models. 
Ethylene (C2H4) fares a bit better: it is produced pre- 
dominantly from C2H3 -I- H2 — C2H4 -I- H (see some 
of the C2H2 production schemes above) and is lost from 
the reverse of this reaction (heading toward C2H2 pro- 
duction), as well as H C2H4 M -> C2H5 + M, to 
eventually form methane. Although the mole fraction of 
C2H4 can exceed that of C2H2 at the C2H2 mole-fraction 
minimum near 0.1-10"^ bar in some of our models, the 
peak C2H4 mole fraction never approaches that of the 
peak C2H2 mole fraction, and the column abundance of 
C2H4 is intermediate between that of C2H2 and C2H6 in 
our HD 189733b and HD 209458b models (see Table [T]). 
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Although our kinetics is far from complete for C3- 
to-Cg species, we do attempt to track the production 
and loss of some complex hydrocarbons. The dom- 
inant C3 hydrocarbon in our models is actually the 
C3H2 radical. This result may be due to our insuf- 
ficient knowledge of CaHj; kinetics, and we encour- 
age further laboratory and theoretical studies of the 
fate of CsH^ species under the low-pressure, high- 
temperature, reducing conditions in the relevant re- 
gions of extrasolar-planet atmosph eres. Both CH and 
C readily insert into C2H2 (e.g., iBerman et al.l [l982l: 
Canosa et al.l 119971: iHaider fc Husain 19931 Kaiser et all 
20011: iCasavecchia et al.l "|2001; Meb el et al.ll2007[ ). and 
these insertion reactions likely dominate the production 
of C3 hydrocarbons on extrasolar giant planets, but the 
fate of the C3H2, C3H, and C3 radicals produced from 
these reactions is unclear. We implicitly assume that the 
C3 and C3H radicals that are produced from the reaction 
of C with C2H2 quickly end up as C3H2 — an assump- 
tion that may only be correct at altitudes where H2 is the 
dominant atmospheric constituent. The main loss pro- 
cess for C3H2 in our models is H C3H2 -I- M — > C3H3 
-|- M. The resulting propargyl radicals can be photolyzed 
to recycle the C3H2, can react with atomic H to reform 
C2H2, can react with H2 or H to form C3H4, or can react 



with other C3H2; radicals to form Cg hydrocarbons and 
eventually benzene. 

Although benzene is not a major constituent in our 
models, it can reach mole fractions of 0.1-1 ppb in the 
few-microbar region in some of our HD 189733b mod- 
els; however, benzene does not survive in the warmer 
HD 209458b models. While ppb levels are by no means 
sufficient for CgHg condensation on HD 189733b, it does 
suggest that polycyclic aromatic hydrocarbon (PAH) for- 
mation could possibly occur at high altitudes on cooler 
transiting exoplanets, which could ultimately lead to the 
for mation of a "soot" a erosol layer, as has been suggested 
bv lZahnle et al.l (|2011h . In contrast, PAH formation will 
be strongly suppressed on planets with strong strato- 
spheric thermal inversions (see also iLiang et al.l [2004) . 
Our neutral chemistry is not particularly conducive for 
soot formation, even under HD 189733b conditions, but 
it is possible that photoionization of N2 or other species 
could lead to the production of high-molecular-weight 
ions and photochemically produced a,erosols through a 
Titan-like ion-chemistry process (e.g.. | Wa ite et al. 20071 
[2OOI llmanaka fc Smith|[2007t iVuittolTeral.. 2007, ,200a) , 
although it is also possible that the high abundance of 
H2O and/or CO and the low abundance of CH4 at high 
altitudes on highly irradiated exoplanets could short- 
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circuit this process. The potential for ion chemistry play- 
ing a role in the formation of complex species in extra- 
solar giant planet atmospheres deserves further study. 

Unlike the case for CO and H2O, methane can be ef- 
ficiently destroyed and not recycled in the upper atmo- 
spheres of HD 189733b and HD 209458b. Some of the 
carbon ends up in C2H2 and other hydrocarbons, as well 
as CO and C, but cross reactions with nitrogen species 
become a key mechanism for removing CH4 from the up- 
per atmosphere. Hydrogen cyanide is the dominant prod- 
uct of this coupled carbon-nitrogen chemistry. From a 
column-integrated standpoint, the dominant photochem- 
ical scheme that permanently converts CH4 into HCN in 
the stratospheres of our exoplanet models is 

2(H20 + ^ U + OU) 
2(0H + H2 ^ H2O +H) 
NH3 + H ^ NH2 + H2 
CH4 + H ^ CH3 + H2 
NH2 + H ^ NH + H2 
NH + H ^ N + H2 
N + CH3 ^ H2CN + H 
H2CN + H ^ HCN +H2 
Net : CH4 -I- NH3 ^ HCN + 3H2. 

(14) 

Hydrogen cyanide then becomes a major disequilib- 
rium product in exoplanet atmospheres and can replace 
methane as the second most abundant carbon-bearing 
species at high altitudes. 

As is the case for oxygen compounds, the photochem- 
istry of carbon compounds is more effective for cooler ex- 
oplanets, both because of the larger abundance of CH4 
from transport-induced quenching and because the lower 
temperatures allow disequilibrium products to survive 
more readily. The dominant carbon-bearing photochem- 
ical products are HCN, C2H2, C, and CH3, with some 
CO2 being produced at high altitudes. 

3.5. Photochemistry of nitrogen species 

Nitrogen photochemistry parallels that of carbon 
chemistry: N2, like CO, is a high-altitude source, 
and NH3, like CH4, is a critical transport-quenched 
species that drives photochemistry throughout the bulk 
of the stratosphere. Figure |6] shows vertical profiles 
of the major nitrogen-bearing species in our nominal 
thermo/photochemical kinetics and transport models. 
At high altitudes, nitrogen photochemistry is initiated 
by N2 photolysis through schemes such as 

CO + hiy^C + 

+ H2 ^ OH + H 
N + OH ^ NO + H 
C + NO ^ CN 
CN + H2 ^ HCN + H 

2 H + M ^ H2 +M 

Net : CO N2 -h H2 ^ HCN + N + + R, 



where the resulting HCN and N are significant high- 
altitude products. On the cooler HD 189733b, additional 
speculative processes involving HC3N rather than NO 
might contribute at high altitudes: 

CO + hi^ + 
N2 + /ii' ^ 2 N 
C -f- C2H2 C3H2 
N + C3H2 ^ HC3N -I- H 
H + HC3N ^ CN + C2H2 

CN H2 HCN + H 

Net : CO + N2 + H2 ^ HCN + N + + H. 



Molecular nitrogen is quite stable on HD 189733b and 
HD 209458b. Although N2 can combine with atomic 
H to form NNH, thermal decomposition of this unsta- 
ble species acts to recycle the N2 throughout the strato- 
sphere. At high altitudes, the strong N=N bond can be 
destroyed by photolysis or attacked by atomic C, H, and 
O to form atomic N, which then can become the domi- 
nant nitrogen species at high altitudes. However, N2 is 
also recycled through schemes such as 

N2 + /li^ ^ 2 N 
N -I- H2 ^ NH + H 
N + NH ^ N2 + H 
Net : H2 ^ 2 H 

and 

N + OH ^ NO + H 
N + NO ^ N2 + 
+ H2 ^ OH + H 
Net : H2 ^ 2H, 

where catalytic destruction of H2 is the net result of both 
recycling schemes. Throughout much of the stratosphere, 
N2 is shielded from photolysis by a larger surrounding 
column of H2 . 

Ammonia has much weaker bonds and is less stable 
than N2, so much of the interesting nitrogen chemistry 
in our models involves NH3. Throughout the bulk of 
the stratosphere, NH3 destruction leads to HC N pr oduc- 
tion through scheme (14) discussed in section |3^ This 
scheme is efficient enough — and HCN is recycled read- 
ily enough — that HCN can take over from NH3 and 
CH4 as important nitrogen and carbon carriers in the 
upper stratosphere. The main loss process for HCN in 
hot- Jupiter stratospheres is the reaction H -I- HCN 
CN + H2, with a much lesser contribution from HCN 
photolysis to form CN -I- H2 . The CN produced by these 
processes reacts with H2 to reform HCN, so HCN recy- 
cling is prevalent. 

Complex nitriles are also formed from photochemistry 
on cooler extrasolar planets like HD 189733b. For exam- 
ple, some of the CN produced in the H -I- HCN reaction 
can react with C2H2 to form HC3N, or some of the CH3 
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produced from hydrogen abstraction of methane can re- 
act with HCN to form CH3CN. These nitriles tend to be 
unstable in a background bath of H, and their destruc- 
tion pathways tend to reproduce the HCN, so they are 
nowhere near as abundant as HCN in our HD 189733b 
models. However, as with benzene and potential PAHs, 
further nitrile chemistry may produce condensible prod- 
ucts, particularly on colder exoplanets. On the other 
hand, complex nitriles will be very unstable on planets 
like HD 209458b that have a strong stratospheric thermal 
inversion. Other carbon-nitrogen bonded species, includ- 
ing imines and amines, are relatively unimportant in our 
models, although the H2CN that is produced during the 
course of scheme (14) can undergo three-body recombi- 
nation with H to form some CH2NH that persists in our 
cooler atmospheric models. Note that nitrogen hydrides 
like hydrazine (N2H4) are unimportant in our models, 
unlike the situation on the much colder Jupiter. 

Cross reactions between nitrogen and oxygen species 
can also occur in exoplanet stratospheres, leading to the 
production of species like NO and HNCO. The dominant 
schemes producing and destroying NO were discussed in 
section l373l HNCO is produced predominantly through 
the reaction NH2 CO ^ HNCO + H, where the NH2 
derives from hydrogen abstraction of ammonia by atomic 



H. The dominant loss for HNCO is reaction with atomic 
H in the reverse of the main production reaction. 

The dominant nitrogen-bearing photochemical prod- 
ucts in our exoplanet models are HCN, N, NO, and NH2. 
As with oxygen and carbon photochemistry, cooler plan- 
ets tend to have more abundant and more varied dise- 
quilibrium nitrogen species than warmer planets. This 
tendency results from the larger abundance of photo- 
chemically active, transport-quenched species like NH3 
on cooler planets and from the greater stability of the 
photochemical products at low temperatures. 

3.6. Sensitivity to temperature 

The atmospheric thermal structure can strongly affect 
the abundance of disequilibrium constituents. As shown 
in Fig. [21 high temperatures in the stratospheres of tran- 
siting hot Jupiters can allow photochemically produced 
species to chemically react with other atmospheric con- 
stituents to drive the composition back toward equilib- 
rium. Disequilibrium compositions are thus difficult to 
maintain on planets like HD 209458b that have strong 
thermal inversions. In contrast, thermochemical kinetics 
is inhibited when stratospheric temperatures are cooler, 
so that photochemically produced species are more likely 
to survive on HD 189733b and other exoplanets without 
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stratospheric thermal inversions. Temperatures in the 
deep atmosphere also affect the abundance of disequilib- 
rium constituents through controlling where transport- 
induced quenching begins to operate. Fig. [7] illustrates 
the importance of this effect for key transport-quenched 
species like CH4, NH3, and HCN. 

In Fig. [71 all models have the same constant-with- 
altitude eddy K^z — 10^ cm^ s~^ profile, but the thermal 
structure is different for each model (see Fig. [2] for the 
adopted temperature profiles). In general, colder tem- 
peratures at depth imply that transport will be able to 
compete with thermochemical kinetics at greater pres- 
sures; the quench point will therefore be deeper on 
colder planets, and the mole fraction of the quenched 
species will be larger. This effect can be seen with the 
HD 189733b models shown in Fig. [71 where the cold- 
est western-terminator-average model has a CH4 quench 
point near the ^7 bar, ~1485 K level, whereas the 
warmest projected-area-weighted dayside-average mod- 
els have quench points near 1.9 bar, ~1522 K. The re- 
sulting CH4 quenched mole fraction is a factor of ~10 
greater for the cooler western-terminator model than for 
the warmer dayside-average model. Ammonia quenches 
at even greater pressures and temperatures, at which 
point the temperature profiles are more uniform (see 
Fig. [5]) , so that the differences between models are less 
pronounced. The HCN profiles are more complicated be- 
cause the quenched abundances of both NH3 and CH4 
affect the HCN abundance, and HCN is photochemi- 
cally produced at high altitudes, but Fig. [71 suggests 
that HCN first quenches in the 10-30 bar regio n, where 
the GCM simulations of iShowman et all ()2009[ ) exhibit 
western-terminator-average temperatures that are actu- 
ally warmer than the dayside-average temperatures at 
these pressures. As such, the western-terminator-average 
model has a smaller quenched HCN abundance than 
the dayside-average models . The 1-D "4pi" profile from 
iFortnev et all (l2006all201Cil) is the coolest of all our mod- 
els at 10-30 bars, resulting in the largest quenched HCN 
mole fraction in this model. Given the lack of constraints 
on temperatures and eddy diffusion coefficients at depth, 
the quantitative results for the abundance of quenched 
species have significant uncertainty in our models; more- 
over, horizontal advection time scales may be shorter 
than chemical time scales at pressures greater than a 
few mbar, acting to homogenize the composition longi- 
tudinally (see section . 

For HD 209458b, the presence of high stratospheric 
temperatures, particularly on the dayside, is the biggest 
factor in controlling the composition, as the disequi- 
librium molecular constituents do not survive in the 
middle and upper stratosphere. The cooler the atmo- 
sphere is at high altitudes, the more likely disequilib- 
rium species are to survive to these altitudes. Note that 
the lack of UV photons in the nightside-average model 
shown in Fig. [71 also helps keep disequilibrium molecules 
present at high altitudes. Even without a strong ther- 
mal inversion at the terminators of HD 209458b, strato- 
spheric temperatures are considerably higher than on 
HD 189733b, enabling thermochemical kinetic processes 
to destroy disequilibrium constituents. Temperatures at 
depth are also much warmer on HD 209458b than on 
HD 189733b, due to the different luminosities of the host 
stars. Transport-induced quenching is therefore less ef- 



fective on HD 209458b than on HD 189733b, and the 
predicted quenched abundances of CH4, NH3, and HCN 
are correspondingly reduced. 

The abundances of HCN and CH4 are particularly sen- 
sitive to temperature, whereas species like H2O, CO, and 
CO2 tend to more closely follow their equilibrium profiles 
such that these species are more sensitive to metallicity 
than to temperature on HD 189733b and HD 209458b. 
This conclusion will differ for colder planets like GJ 436b 
that are in the CH4 rather than CO stability field. 

3.7. Sensitivity to eddy diffusion coefficients 

The model results are also sensitive to the assumed 
eddy diffusion coefficients, again because of the potential 
importance of transport-induced quenching on extraso- 
lar giant planets. Fig. [SI demonstrates this sensitivity. 
For larger eddy diffusion coefficients, transport-induced 
quenching becomes effective at deeper pressures, lead- 
ing to larger mole fractions of quenched species. Larger 
Kzz values also allow molecular constituents to be carried 
to higher altitudes. Both these processes contribute to 
enhanced abundances of photochemically derived prod- 
ucts when atmospheric mixing is more vigorous. For ex- 
ample, when quenched species like CH4 and NH3 are 
more abundant and carried to higher altitudes, photo- 
chemical products like HCN and other nitriles, and C2H2 
and other hydrocarbons can build up over a larger col- 
umn. The H2-to-H conversion level also changes with 
Kzz, which can have a significant effect on all molecular 
abundances. Minor heavy constituents like complex hy- 
drocarbons and nitriles are particularly affected by Kzz, 
as these species often depend nonlinearly on the abun- 
dance of parent molecules like CH4 and NH3 . Note from 
Fig. m that if K^z is sufficiently high, there is even a 
quench point for H2O and CO that differs from the equi- 
librium solution. 

General ci rculation models for highly irradiated hot 
Jupiters fsee lShowman et aLll2009l fShowman et al.ll20Tol 
and references therein) suggest relatively strong verti- 
cal winds, which according to our assumptions could 
lead to potentially vigorous mixing at pressures less than 
~100 bar [i.e., Fig. [ij. We have assumed that con- 
vection in the adiabatic portion of the atmosphere at 
grea ter depth als o contributes to vigorous vertical mixing 
(see IStond Il976[) . Transport-induced quenching is then 
very important, and quenched disequilibrium species and 
their photochemical products are abundant. In con- 
trast, Youdin & Mitchell (2010) suggest that weak tur- 
bulent mixing in extrasolar giant planet atmospheres 
could slow the cooling rate of the atmosphere and might 
explain the l arge, "inflated" rad i i of s ome hot-Jupiter 
atmospheres. lYoudin fc Mitchellj ()2010f) further suggest 
that Kzz values greater than 10^ cm^ s~^ would over- 
inflate hot- Jupiter radii beyond what can be supported 
by observations — they favor much smaller Kzz val- 
ues of (1-5) X 10^ cm-^ s^^ at ^550 bar, or in general 
Kzz < 10^ cm'^ s"^ at P > 100 mbar for tempera- 
tures of Tdeep = 1500 K in the deep isothermal radia- 
tive region above the optically thick region of the at- 
mosphere, similar to the few-bar-to-few-hundred-bar re- 
gio n for our HD 1897 3 3b te mperature proflles. With 
the lYoudin fc Mitchellj (|201 J ) formalism, Kzz scales as 
as Tl so that Kzz on HD 209458b could be slightly 
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A color version of this figure is available in the online journal. 
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higher by a facto r of 2-3. If eddy d if fusion coefficients 
were as small as lYoudin fc Mitchelll (|20Tol) suggest in 
the ~2-500 bar region, then transport-induced quench- 
ing would be significantly suppressed in our HD 189733b 
and HD 209458b models, and the quenched abundances 
of CH4, NH3, and HCN would be quite low — more con- 
sistent with what is predicted from thermochemical equi- 
librium — and these species (including methane) would 
likely not be observable. Although waves propagating 
up from the deep interior and down from the weather 
layer may keep K^z from being as low as is suggested 
bv lYoudin fc Mitchel 42010), the idea of a "mechanical 
greenhouse" is an interesting one. Spectral observations 
that can unambiguously confirm the presence of CH4, 
NH3 , or HCN at levels in excess of equilibrium could help 
constrain the strength of turbulent mixing in extrasolar- 
giant-planet atmospheres. 

Given that the NH3 quench point is likely in the adi- 
abatic region of the deep atmosphere, it is interesting to 
speculate that the NH3 abundance might eventually help 
provide a measure of the entropy of the adiabat, which 
is currently known only crudely from the planet's radius. 
For a given mass, planets with larger radii have higher- 
entropy adiabats and lower-pressure radiative-convective 
boundaries. Although the entropy of the adiabat will 
certainly affect the NH3 quench point, the overall uncer- 
tainties in the eddy diffusion coefficients, in the planetary 
metallicity, and in the kinetics of NH3 quenching might 
make such a suggestion impractical. 

3.8. Photochemistry and transport time constants and 
their implications 

Because our models remain at a single solar zenith an- 
gle for the long duration of the simulation, they are not 
true analogs for the situation on a real planet, even if 
the planet is tidally locked and has one hemisphere al- 
ways facing the star. Real planets have meteorology, and 
dynamical models indicate that strong winds will be act- 
ing on close-in transiting exoplanets to transport s pecies 
zonal ly, meridionally, and vertically (see Showman et al.l 
[201^ for a review of exoplanet dynamical models); zonal 
jets are predicted to be particularly strong, such that ma- 
terial is transported longitudinally across a hemisphere 
with a horizontal t ransport time scale Tdyn /i ~ 2 x 10^ 
s on HD 209458b (|Cooper fc Showma n 2006). A parcel 
of gas will therefore experience variable external forc- 
ing from the continuously changing irradiation angles, 
and those effects can influence the chemistry. A full 3- 
D model of the radiative, dynamical, and chemical cou- 
pling needed to describe this situation would be com- 
putationally prohibitively expensive, but simplifications 
such have been described by Cooper & Showman (2006) 
can be made to track the key chemical effects. Although 
we do not attempt these kinds of simulations with our 
simple 1-D models, we can at least examine the time 
constants in the system and make some comments about 
what the time-constant analysis implies about the com- 
position of mo re realis ti c plan etary analogs. 

Cooper fc ShowmanI (l2006[) convincingly argue that 
the large CO abundance that is predicted for day side 
HD 209458b will not be chemically converted to CH4 
at the terminators or in the nightside stratosphere, de- 
spite the fact that CH4 would be favored in chemical 
equilibrium. Vertical transport time scales are simply 



much shorter than the chemical conversion time scale be- 
tween CO and CH4. In addition, note that just as there 
is transport-induced quenching in the vertical direction, 
horizontal transport-induce d quenching can also operate 
(iCooper fc Showman"2006') . Will this horizontal quench- 
ing affect photochemically produced species? One might 
particularly want to know whether HCN will reconvert 
into CH4 and NH3, whether C2H2 will reconvert into 
CH4, and whether H will reconvert to H2 at high alti- 
tudes on the nightside of HD 189733b and HD 209458b, 
and whether vertical transport-induced quenching will 
lead to different CH4, NH3, and HCN mole fractions on 
the day side as compared with the terminators (as pre- 
dicted by our models) or whether horizontal transport- 
induced quenching will homogenize the quenched mole 
fractions in the stratosphere. 

In general, the horizontal advection time scale over a 
planetary radius is comparable to or less than the verti- 
cal advection time scale over a scale height, suggesting 
that if temperatures are low enough in the upper tro- 
posphere and lower stratosphere to shut down CH4 o 
CO, NH3 o N2, and CH4 O NH3 o HCN intercon- 
version, then horizontal quenching will also be in op- 
eration in these regions. As such, zonal winds would be 
expected to operate efficiently enough to homogenize the 
CO, CH4, H2O, N2, and NH3 abundances in longitude, 
and there should be no terminator- vs. -dayside differences 
in the mole fractions of these species in regions where 
those interconversion reactions dominate. However, in 
the middle and upper stratosphere, the disequilibrium 
that results from the absorption of UV photons has a 
different set of quench reactions operating under much 
different conditions. Photochemically produced species 
tend to have shorter lifetimes, and vertical transport 
time scales are also expected to be shorter at high al- 
titudes. Zonal wind s may also increase w ith altitude in 
low- latitude regions (jShowman et al.|[2009D . but horizon- 
tal transport time scales eventually become longer than 
the photochemical time constants. Therefore, longitudi- 
nal differences are more likely to be maintained at higher 
altitudes (i.e., pressures less than '-^O. 01-0. 001 mbar, de- 
pending on the species and model). 

In our nominal dayside HD 189733b atmosphere for 
instance (see Fig. [5]), the HCN that replaces CH4 and 
NH3 at pressures less than a few mbar is produced 
through scheme (14). This scheme does not operate ef- 
fectively at pressures greater than a few mbar, and ver- 
tical transport dominates the behavior of the HCN pro- 
file at these pressures such that the HCN that was pro- 
duced at higher altitudes diffuses with a constant flux 
to the lower atmosphere; however, the HCN photochem- 
ical lifetime is only ~2 hr at 10~^ mbar and less than 3 
min at 10"'^ mbar, compared to a zonal transport time 
scale that is likely gre ater than 20 hr at these pressures 
(jShowman et al.ll2009[ ). The photochemical lifetime of 
C2H2 is even shorter, suggesting that longitudinal vari- 
ations in HCN, CH4, NH3, and C2H2 could be main- 
tained at high altitudes on HD 189733b, such that CH4 
and NH3 could recover at high altitudes at night. The 
total column abundances of CH4 and NH3 would then be 
larger on the nightside than on the dayside, particularly 
at pressures less than '^l mbar. 

It is interesting to speculate what might happen on 
HD 189733b when the high-altitude CH4 that comes 
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back at night is transported by the zonal winds to 
the western (dawn) terminator (assuming atmospheric 
super-rotation, as is expected from the dynamical mod- 
els). As the gas flows into the starlight again, the high- 
altitude CH4 will strongly absorb in the ultraviolet and 
at certain visible and near-IR wavelengths where the 
methane bands have some optically thick lines. Does 
this expected "pulse" of high-altitude absorption on the 
dawn limb affect the high-altitude thermal structure or 
the CH4 fluorescent emission behavior enough to cause 
emission in these optically thick lines? If the heating 
rates are greater than the chemical destruction rates, 
if there is sufficient methane at high altitudes, and if 
there is a lag in the cooling rates, then such a scenario 
might help explain the ap parent 3.25 /im emission seen 
in the iSwain et al.l ()2010t ) ground-based observations of 
HD 189733b. We note that if this process were really 
operating, the emission in the vs methane band should 
be more obvious in the transit of the dawn terminator 
than the dusk terminator because more methane will be 
present at high altitudes on the dawn terminator (again, 
assuming atmospheric super- rotation). Even if cooling 
rates were rapid, the sudden pulse of photons at high alti- 
tudes would provide energy needed to drive the chemical 
reactions or fluorescence and might excite any molecules 
present. 

This process would work on HD 209458b as well, but 
it would be less apparent due to the overall smaller 
quenched abundance of CH4 and the correspondingly 
smaller column abundance of CH4 at high altitudes. On 
planets like GJ 436b that are in the CH4 stability regime 
at the deep CO o CH4 quench point, the effect could 
be particularly important. Because CH4 is likely rela- 
tively abundant at high altitudes on both terminators on 
GJ 436b, we suggest that high-spectral-resolution mod- 
eling of radiative absorption and emission on GJ 436b 
might be called for to determine whether portions of the 
V3 band would be expected in emission rather than ab- 
sorption on this planet (due to a high-altitude thermal 
inversion from UV, visible, and near-IR absorption from 
CH4), potentially explaining the apparent lack of CH4 
absorption at 3.6 um in eclipse observations of this planet 
([Stevenson et al.lF2010[ ). These ideas are quite specula- 
tive at this point but deserve further quantitative study. 
The possibility of dawn-dusk differences in the transit 
spectra could potentially be expl ored with the James 
Webb Space Telescope (JWST; see lFortnev et al. 2010). 

4. IMPLICATIONS WITH RESPECT TO OBSERVATIONS 

Figure ini illustrates how the predicted HD 189733b and 
HD 209458b spectra from the primary transits are altered 
by the presence of the disequilibrium products derived 
in our models. Different radiative-transfer models in the 
literature are producing different predicted transmission 
spectra for exoplane t s, for reasons that ar e unclear (cf. 
iShabram et all |2QTiI: iBeaulieu et aI|[20Tll ). For Fig. M 
the synthetic transmission s pectra are calculated as is 
described in iFortnev et al.l (|2010l ) and ISh abram et al.l 
(poll.) for the terminator-average temperature profiles 
shown in Fig. [2] The calculations assume either a 
thermochemical-equilibrium composition or use the re- 
sults from our thermochemical and photochemical kinet- 
ics and transport models with the nominal K^^ profiles. 
The approximate pressure levels from which the spectral 



features originate are discussed further in IFortnev et al.l 
(|2008aD . For HD 189733b, Fig. |9] shows clear excess ab- 
sorption from the disequilibrium model at wavelengths 
near ~ 2. 1-2. 5 fim and ^2.9-4 fim due to CH4 and to a 
lesser extent due to HCN near 3 /im, additional absorp- 
tion in the ~7-9 fim region due largely to CH4 (with 
HCN contributing near 7 /xm), and additional absorp- 
tion in the ~9-15 /xm region due to NH3 (centered near 
10.5 /im) and HCN (centered near 14 fim), with a much 
lesser contribution from C2H2 (centered near 13.6 /im). 
Note that the column abundance of C2H2 in our nomi- 
nal model is small enough that acetylene absorption fea- 
tures are not very prominent. For HD 209458b, spectral 
differences between the equilibrium and disequilibrium 
models are much more muted, as the overall composi- 
tional changes are slight between the equilibrium and 
disequilibrium models at the pressure levels at which 
the observations are sensitive. Our modeling suggests 
that photochemistry and transport-induced quenching 
will have only a minor effect on the spectral properties 
of the hottest "hot Jupiters" or those with strong strato- 
spheric thermal inversions, as the high temperatures en- 
able rapid kinetics that can drive the atmosphere back 
toward equilibrium. 

For HD 189733b, however, the disequilibrium chem- 
istry resulting from photochemistry and transport- 
induced quenching measurably alters the composition 
such that the disequilibrium model shows additional 
absorption and a greater wavelength-dependent con- 
trast overall than the equilibrium model. Although 
the near-infrared absorption band positions predicted 
by the model are generally consistent with observa- 
tions, both the equilibrium and disequilibrium models 
apparently underestimate the near-IR band depths and 
the overall amplitude of the wavelength-dependent vari- 
ations — a fact previously d iscussed for equilibrium 
models by IFortnev et al.l ()2010l ) . Even when disequilib- 
rium constituents are consider ed, the model fit to the 
near-IR ggT/NI CMOS data of' Swain et all ()2008bD and 
iSing et all (|2009l ) is not significantly improved, as wa- 
ter absorption and H2 collision-induced absorption dom- 
inate the spectral behavior at near-IR wavelengths, and 
these opacity sources remain roughly the same in the 
two models. The main exceptions to this statement are 
wavelength regions centered near 2.3 /im and 3.3 /im, 
where methane absorption can have a significant effect. 
The additional methane in the disequilibrium models 
for HD 18 9733b improves the fit in the 2.3 um. region 
(see also ISwain et al.l l2008bt iMadhusudhan fc Seageii 
I2009D but degrades the fit to the IRAC 3.6 /xm data 
(e.g.. 'Bea ulieu et al.|[2008t iDesert et~aIll2009D . At mid- 
IR wavelengths, the models also appear to underesti- 
mate the spectral contrast between different wavelengths, 
and neither model does a good job of reproducing the 
Spitzer /IRAC photometric data, although the data have 
large error bars, and different groups using different anal- 
ysis procedures for the same datasets do not always agree 
on the observe d fiux in the IR AC bands. As was dis- 
cussed by For tnev et al.l (|2010l ). higher metallicity mod- 
els with their increased CO and CO2 abundances would 
improve the fit to the relative absorption strengths of 
the 3.6, 4.5, and 5.8 ii u IRAC band data points from 
the IDesert et al.l (|2009| ) analysis, but other model-data 
comparison problems would still remain. Note that we 
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Fig. 9. — Synthetic transit spectra calculated for our HD 189733b (left) and HD 209458b (right) thermochemical and photochemical 
kinetics and transport models (red curves) that assume a terminator-average thermal structure and nominal K^^ profile, compared with 
synthetic transit spectra from the thermochemical equilibrium models (black curves) for the same assumed thermal structure. All models 
assume a Ix solar composition. Absorption depth is calculated as the square of apparent planct-to-star radius ratio. Observations for 
HD 189733b are shown as data points with as sociated error bars: green, Swain et al. (2008b) //ST/NICMOS; purple. Sing et al. (200S|) 
FST/NICMOS; light blue. Desert et al.l l(2009h Spitzer /IKAC; darker blue, Beaulieu et al. (2008) for Spitzer/IKAC 3.6 ^lm and 5.8 ^tm, 
[Knutson et al. (2007) for Spitzer /mXC 8 /xm, and Knutson et al. (2009) for Spitzer /mAC 247tni; orange, Agol et al. (2010) Spitzer /m,KC 
8 /xm. For HD 209458b, the green data points represent the Spitzer AC data of . Beaulieu et al.. (.2009. ) and the blue data point at 24 /^m 
represents the average of the 5pi<zer/MIPS values from [Richardson et al.l (|2006l ) and H. Knutson (2009, private communication). A color 
version of this figure is available in the online journal. 
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Fig. 10. — Synthetic transit spectra (apparent planetary radius in units of Jupiter radii) calculated for HD 189733b for different as- 
sumptions about the thermal structure (left) and the eddy diffusion coefficient profile (right). In the left panel, the gray curve represents 
a model with the west-terminator-average temperature profile and an equilibrium composition, the red curve represents a model with 
the west-terminator-average temperature profile and a composition indicated by our disequilibrium model with Kzz = 10^ cm'^ s~^, the 
black curve represents a model with the east-terminator-average temperature profile and an equilibrium composition, and the orange curve 
represents a model with the east-terminator-average temperature profile and a composition indicated by our disequilibrium model with 
Kzz = 10® cm^ s~^. In the right panel, the models assume a dayside-average thermal structure: the red curve is the equilibrium model, 
and the other curves represent disequilibrium models with assumed constant K^^ values of 10^ (green), 10** (dark blue), 10® (cyan), 10^" 
(magenta), and 10^^ (orange) cm'^ s~^. A color version of this figure is available in the online journal. 



have not recalculated the thermal profile that will result 
from the additional disequilibrium opacity sources, and 
the expected altered thermal structure may also affect 
the observed spectrum. 

As an example of the sensitivity of transmission spectra 
to the thermal profile, Fig. [10] shows how the synthetic 
spectra for HD 189733b change for different assumptions 
about the thermal structure at the terminators — par- 
ticularly the western vs. eastern terminator-average pro- 
files from Fig. [2] Fig. [TU] also shows the sensitivity of the 
transmission spectra to different assumptions about the 
eddy diffusion coefficient profile, for assumed constant 
Kzz profiles ranging from 10^ to 10" cm s . The 



thermal structure affects such things as the atmospheric 
scale height, which describes how extended the atmo- 
sphere becomes, so that a colder planet will have less of 
an apparent cross section during transit for all other fac- 
tors being equal. As such, the colder western-terminator 
model exhibits a smaller apparent radius than the eastern 
terminator model. Note that the spectral changes caused 
by the equilibrium vs. disequilibrium models have less 
of an overall impact on the apparent absorption during 
transit than assumptions about the thermal structure. 
The eddy diffusion coefficient profile affects the extent 
to which heavy molecular constituents can be mixed to 
high altitudes on an H2-dominated planet; in general, the 
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larger the eddy diffusion coefficient, the higher the alti- 
tude to which the molecular constituents can be carried, 
which allows photochemical species to be produced over 
a larger column of the atmosphere, leading to larger col- 
umn abundances and increased absorption by the photo- 
chemical products. Eddy diffusion coefficients also affect 
the abundances of the species controlled by transport- 
induced quenching, as is discussed in sections 13.21 and 
[3Jl 

Synthetic emission spectra of HD 189733b and HD 
209458b for dayside-average, secondary-eclipse condi- 
tions are shown in Fig. [TTJ For these calculations, lo- 
cal thermodynamic equilibrium (LTE) conditions are as- 
sumed, and the effects of the important spectrally active 
gases H2O, CO, CH4, CO2, HCN, NH3, C2H2, and C2H6 
are considered. The opacities of these gases are garnered 
from different sources. The hot water line parameters 
from the latest November 2010 version of the HITRAN 
line list (iRothman et al.l I2009D are a dopte d, using the 
original data of iBarber et al.l ()2006[ ) and iZobov et al.l 
([2008). Absorption by CO2 is calculated fro m the HI- 
TRAN hot CO2 line hst (Tashkun et aL||2003D, as is that 
of CO and NH3 (Jlothnian et al. 200^" 

and the methane 

li ne parameters derive fro m the hot CH4 line parameters 
of iFreedman et all (|2008[ ). Note that for both the trans- 
mission and emission modeling, the opacities of HCN and 
C2H2 — and really all species other than H2O, CO2, and 
CO — are likely underestimated because the line parame- 
ters were derived from cool, low-temperature conditions. 
Absorption coefficients are calculated using a line-by-line 
technique every 0.004 cm^^ wavenumbers to provide the 
k-coefficients that are included in the radiative-transfer 
calculations. For the calculations of the emission spectra, 
the atmosphere is divided into 80 vertical grid points that 
extend from 10^^ to 10 bar. The effects of particulates 
are excluded. The calculations explore the spectroscopic 
effects of the compositions that are indicated by the dif- 
ferent chemical models, while the temperature structure 
is held fixed. The depths at which the various emission 
signatures originate can be seen from the contribution 
functions shown in Fig. 1121 

The synthetic emission spectra shown in Fig. [TT] il- 
lustrate that spectral signatures from the disequilibrium 
model are noticeably altered from the equilibrium predic- 
tions for the cooler HD 189733b but not for the warmer 
HD 209458b. The predicted eclipse depths are only 
slightly modified at long wavelengths due to the dise- 
quilibrium chemistry on HD 209458b. That conclusion 
changes for HD 189733b, as the the key disequilibrium 
constituents CH4, HCN, and NH3 are more abundant 
and have a larger impact on the predicted spectrum. On 
HD 189733b, the disequilibrium-model spectra diverge 
from the equilibrium spectra in the 7-9 fim range due to 
the z/4 CH4 band at 7.7 /xm, in the ~3-4 /xm and ^2.1-2.5 
/Lim ranges due to the 3.3 and 2.3 ^m bands of CH4, in 
the 5.5-6.6 /xm range due to the 1/4 band of NH3 at 6.15 
/zm, in the 8-12 /im range due to the 1^2 band of NH3 at 
10.5 /im, and in the 13-15 /im range due to the 1^2 band 
of HCN at 14 /im. The individual contributions from the 
different disequilibrium species are better illustrated in 
Fig.Hl 

Our emission spectra for both HD 189733b and 
HD 209458b do not compare well with the observed 
fluxes during secondary eclipse (e.g., iDeming et al.l 
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at 4.5 and 5.8 fim (see, e.g., Knu tson et al.l '20081 are 
not present in the model spectra, which could be due 
to an inaccurate temperature profile adopted for our 
model (e.g., a thermal inversion beginning at higher 
altitudes in our model than on the real planet), or some 
other inappropriate assumptions for the spectral or 
chemical models. For HD 189733b, the model eclipse 
depths are much smaller than has been observed at all 
wavelengths, perhaps resulting from inaccurate model 
assumptions such as the atmospheric thermal structure 
of the planet, the planet's metallicity and radius, or 
the spectrum of the host star (note th at we have used 
the models of iCastelli fc Kuruc3 120041 to simulate the 
spectrum of the host stars, with other stellar an d 
planetary parameters as described in lTorres et al.ll2008() . 
In particular, the model-data comparisons shown in 
Fig. [11] suggest that our assumed dayside atmospheric 
temperatures may be too cold in the 0.01-1 bar region. 
We therefore also calculat e synthetic sp ectra for the 
warmer 1-D "27r" model of iFortnev et al.l (|2006aL IMol l 
(see Fig. [5]); the results are shown in Fig. [TJl along with 
the corresponding contribution functions. Note that the 
magnitude of the eclipse depths are better predicted and 
the disequilibrium effects are reduced with this warmer 
model. 

Our emission spectra for both HD 189733b and 
HD 209458b do not compare well with the observed 
wavelength-dependent fluxes. Rather than dwelling on 
possible reasons, we have focused on the expected differ- 
ences in the spectral signatures of equilibrium and dis- 
equilibrium compositions in this section. Detailed com- 
parisons of our predicted mole fractions with inferences 
from published analyses of transit and eclipse observa- 
tions are deferred to sections [5] and [T] 

5. COMPARISONS WITH OTHER PHOTOCHEMICAL 
MODELS 

We now compare o ur results wi t h tho se from other 
photochemical models. iLiang et al.l ()2003f) were the first 
to develop a photochemical model for hot- Jupiter exo- 
planets. Their main conclusion from their HD 209458b 
modeling was that atomic H is produced in large quan- 
tities from photolysis of H2O and from the subsequent 
reaction of OH with H2 (see our scheme (10) in section 
13.31) . We agree with this general result, although some 
quantitative differences exist between our models. For 
instance, atomic H never replaces H2 a s the dominant 
constituent at high altitudes in the Liang "ei~all (f2003h 
models, for reasons that are unclear. O ur HD 209458b 
mo dels are more comparab le to the ! Garcia Munod ()2007l ) 
and lKoskinen et all ()2010f) models with regard to this H- 
abundance issue. 

O ther quantit at ive d ifferences are apparent between 
the ILi ang et al. (2 003[ ) models and ours. Because 
ILiang et al. (2003) do not consider thermochemical ki- 
netics and do not fully reverse their reactions, they take 
a reasonable approach in beginning their model at a 
base pressure of ^ 1 bar, where equili brium effects will 
be less prevalent. ILiang et al.l (|2003 l) do not consider 
transport- induced quenching of CH4, and their assumed 
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Fig. 11. — Synthetic emission spectra for HD 189733b (left) and HD 209458b (right) for our thermochemical and photochemical kinetics 
and transport models (red curves) that assume a dayside-average thermal structure and corresponding nominal K^^ profile for the planet 
in question (see Figs, [l] and |2]l, compared with synthetic spectra from the thermochemical equilibrium model (blue curves) for the same 
assumed thermal structure. All models assume a Ix solar elemental composition. Spectra are calculated as the ratio of the flux of the 
planet to the flux of the star. Also shown are data points from various secondary-eclipse observations, as labeled. A color version of this 
figure is available in the online journal. 




Fig. 12. — Contribution functions, which give the relative contribution of a given pressure level to the observed emission at any particular 
wavelength, are shown alongside the thermal structure from our dayside-average chemical models for HD 189733b (left) and HD 209458b 
(right). A color version of this figure is available in the online journal. 



1-bar methane abundance is quite a bit lower than ours; 
however, even ignoring this difference, the shape of their 
methane profiles throughout the atmosphere is quite dif- 
ferent from our HD 209458b models. In the lLiang erall 
(|2003[ ) models, the CH4 abundance increases from 1 bar 
to 0.1 mbar because of production of CH4 from CO pho- 
tolysis, whereas our CH4 mole-fraction vertical slopes 
above the quench point are more constant or even neg- 
ative before methane becomes sharply depleted when 
schemes (14) and (12) become effective, indicating a net 
loss rather than net production of methane in the 1~ 
10~^-bar region o f our models . How ever, both our mod- 
els and those of iLiang et"al] (|2003f) suggest that CH4 
does not survive at high altitudes on HD 209458b. In 
our models, CH4 is destroyed by reaction with atomic 
H (see schemes (12) and (14), for example) and the car- 
bon ends up in sp ecies like HCN and C2H2. In contrast, 
ILiang et al.l (|2003f ) state that CH4 photolysis rather than 
reaction with atomic H causes the high-altitude deple- 



tion, but it is not clear from their discussion where the 
carbon ends up. 

Some of these issue s are clarified in the follow-up model 
of ILiang et al.l (|2004f ). where it seems that hydrocarbons 
like C2H2 may replace CH4 as the second most abun- 
dant carbon bearing species (behind CO) at high alti- 
tudes in their HD 209458b models. The ILiang et all 
(|2i00l "Model A" is closest to our HD 209458b west- 
terminator-average model in terms of the CH4 abundance 
at 0.1 mbar, so these two models can be compared in de- 
tail. Our model has some similar general behavior, such 
as the fact that C2H2 is the dominant C2HJ; hydrocar- 
bon from a column-integrated standpoint and that C2H2 
can survive to relatively high altitudes, as compared with 
C2Hg. We also agree that hydrogenation is the major loss 
mechanism for C2Ha; hydrocarbons that prevents com- 
plex hydrocarbons from becoming major constituents on 
HD 209458b. However, there are also some large differ- 
ences between our models in terms of the overall abun- 
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Fig. 13. — S ynth etic emission spectra for HD 189733b, as de- 
scribed in Fig. 111! except the green curve shows how the results 
change when the contribution of NH3 is removed from the kinetics 
and transport model spectra, and the orange curve shows how the 
results change when the contributions of both NH3 and HCN are 
removed from the kinetics and transport model spectra. Note that 
methane opacity is the primary source of the divergence of the or- 
ange curve from the the blue equilibrium curve. The dashed lines 
are blackbody curves for particular temperatures; top curve is for 
1400 K, middle curve is for 1200 K, bottom curve is for 1000 K. A 
color version of this figure is available in the online journal. 

dances and vertical profiles of C2Ha; hydrocarbons. For 
example, our C2H2 profiles tend to peak at high altitudes 
and fall off sharply with both increasing and decreasing 
altitude, whereas their C2H2 profiles are roughly con- 
stant above ~0.1 mbar. In our models, three-body re- 
actions with atomic H that act to destroy C2H2 become 
increasingly effective with increasing pressure, and C2H2 
ends up being more efficiently converted to CH4 in the 1- 
10^'^-mbar region (see scheme (13) in sect ion than 
is apparently the case in the iLiang etlUI ([2004) mod- 
els. Also, C2H2 cannot survive the large H abundance 
above 10 nbar in our model, whereas C2H2 read ily sur- 
vives at those pressures in the ILiang et al.l (200j) model, 
suggesting that H is much less abundant or that C2H2 
conversion into atomic C is m uch less effective in this re- 
gion of the ILiang et al.l (|2004( ) models. In general, C2lla; 
hydrocarbons are much less abundant in the l-10~^-bar 
region of our HD 209458b models, suggesting that our 
kinetics mo re efficiently conve rts these species to CH4 
than in the ILiang et al.l ()2004l ) models (see scheme (13) 
in sec tion 13.41). 

The lLiang et all (|2001 . |2004D models were updated sig- 
nificantlv bv iLine et al. ()2010[ ). who consider methane 



quenching for the first time and expand th e range of 
high-t emperature reactions in their kinetics. iLine et aP 
(I2OIOD use the time- scale arguments developed by 
iPrinn fc Barshavl (|1977D (see section to estimate the 
quenched CII4 mole fraction on HD 189733b. Their 
resulting quenched mole fraction for CH4 — 4 x 10~^ 
— is reasonably consistent with the results of our ki- 
netics/transport model (see Fig. [5]), but their reason- 
ing in deriving this value is marred by the fact that 
they incorrectly reverse their assumed ra te-li i niting CO 
CH4 conversion step (see [Bezard e t al.l 120021 and 



I4 conversion s tep (see 
iVisscher et alll2010bt iVisscher fc Mosesi 120111 for infor- 
mation on correct reversal procedures for this assumed 



rate-limiting reaction) and by the fact that the rever- 
sal was not necessary in the first place since it is CH4 
CO quenching that matter s for HD 189733b . Given 
these issues and the fact that ILine et all (j2010( ) have a 
different rate-limiting step for CH4 quenching than we 
have derived in our m o del, a s well as the fact that they 
do not use the iSmithI (|1998[ ) procedure when calculat- 
ing the transport time scale, it is remarkable that their 
quenched mole fraction for CH4 is so close to that in our 
models. 

In any case, the procedure used by ILine et all ()2010l ) 
once the CH4 quench point has been established is per- 
fectly valid — it may even be preferable to full ther- 
mochemical and photochemical kinetics and transport 
models in some instances. As an example, the quench 
kinetics for CH4 has some uncertainty associated with it 
due to uncertainties in the kinetic rate coefficients, the 
eddy diffusion coefficients, and the atmospheric thermal 
structure. If the CH4 abundance for a transiting exo- 
planet were unambiguously determined from transit or 
eclipse observations, it might be better to take this in- 
formation into account and begin photochemical calcula- 
tions at some point above the quench level (i.e., put the 
bas e of the model ju st above the CH4 quench level, as in 
the lLine et al.|[20lol models). This procedure would work 
best for C-H-0 chemistry investigations, as nitrogen- 
bearing and sulfur-bearing constituents have yet to be 
detected in hot- Jupiter atmospheres. 

The quenched CH4 mole fraction in the iLine et aD 
f|201Qi) daytime model is most similar to our terminator- 
average HD 189733b model (see Figs. H & [5]), and the 
two can be directly compared. In contrast to equilibrium 
predictions, both of our models show a greatly enhanced 
H abundance in the middle and upper stratosphere due 
to scheme (10), both show CH4 depletions at high alti- 
tudes with corresponding significant C2H2 enhancements 
at these altitudes, both show small enhancements in the 
CO2 mole fraction at high altitudes due to CO and H2O 
photochemistry, and both show that photochemical pro- 
cesses have little net effect on the H2O and CO profiles. 
However, there are also some significant differences be- 
tween our rn o dels that ar e worth mentioning. As with the 
ILiang et all H2003I . I2004D photochemica l models, the H 
production rate in the lLine et al.l ()2010f ) models is large, 
but H never exceeds the H2 abundance at pressures < 1 
nbar, as it does in our HD 189733b models. Again, the 
reason for this difference between our models is unclear, 
but it seems to have something to do with the efficiency 
of scheme (10) and its reverse or with other H2 destruc- 
tion and recycling mechanisms. The fact that we have 
considerably more O an d OH at h igh altitudes in our 
model as compared with ILine et aD (2010) is consistent 
with our greater depletion of H2 and is consistent with 
this suggested cause. Our models also show less net pro- 
duction of CO2 fr om photochemical processes than in the 
ILine et all ()2010D model, as thermochemical processes in 
our model drive the CO2 abundance back toward equi- 
librium. 

Another difference between our model and that of 
ILine et al.l (|2010( ) is the more abrupt and more severe 
depletion in the CH4 mole fraction above ~1 /zbar in our 
HD 189733b model as compared with their model. This 
difference is likely due to coupled carbon-nitrogen chem- 
istry (see scheme (14), for example), as HCN takes over 
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Fig. 14.— Synthetic emission spectra for HD 189733b (left) for the 1-D "27r" thermal structure of ll''ortnev et al.l I l2006al . [20lOt) (i.e., 
the temperature structure is from a 1-D radiative-convection equilibrium model where it is assumed that the incident stellar energy is 
distributed over the dayside hemisphere only), along with the contribution functions for this model (right). The model shown by the red 
curve uses results from our thermochemical and photochemical kinetics and transport model with the nominal Kzz profile, whereas the 
model shown by the blue curve assumes a thermochemical-equilibrium composition. Both models assume a Ix solar elemental abundance. 
Also shown are data points from various secondary-eclipse observations, as labeled. A color version of this figure is available in the online 
journal. 



from CH4 at high ahitudes in our mo del, whereas nitro- 
gen chemistry is not considered in the iLine et all (|20100 
models. A more minor difference between our model and 
theirs is the dominant schem e to form C ^ ^ifr h ydrocar- 
bo ns in the stratosphere. In ILine et aD (|2010( ) and in 
the lLiang et al.l[2003L 120041 models on which these mod- 
els were based, the dominant reaction that produces C2 
species from CH^; species is apparently '^CH2 + '^CH2 
— ?• C2H2 -f 2H, whereas we find the reaction ■^CH2 + 
CH3 — )■ C2H4 + H to be more important from a column- 
integrated stand point. Qua l itative ly, however, the main 
behavior in the ILine et aD (|2010l ) models is similar to 
ours, other than some differences that results from our 
inclusion of nitrogen photochemistry. 

We can also compare our models to the thermochemi- 
cal and phot o chem ical kinetics and transpor t mode ls of 
iZahnle et all (|201lD . Note that lZahnle et all ()2009D dis- 
cuss sulfur photochemistry only, which we do n ot con- 
sider, so we do not compare our results with Zahn le et al.l 
(pool . Direct comparisons between our models and 
those of IZahnle et al.l (|2011[ ) are complicated by the fact 
that they show results only for models with isothermal 
atmo spheres. Compariso ns are also hampered by the fact 
that IZahnle et aD (|20lH ) are working on a revision of 
their model as we are preparing this paper here, so that 
we can only compare with their "version 2" in the arXiv 
archive, rather than with a published version of the mod- 
els. Although we cannot do anything about the second 
problem, we can resolve the first problem — which is a 
significant one — by developing an isothermal model of 
our own to better compare with the Zahnle et al. (2011) 
results. The isothermal model that we develop is for a 
Ix solar composition atmosphere, with T = 1000 K, and 
Kzz at a constant 10® cm^ s~^, for a planet located 0.1 
AU from a Sun-like star, at a fixed solar zenith angle of 
30° (see Fig. [T5|) . This model can be directly compared 
with the results show n in the top ri ght figure of Fig. 4 of 
the arXiv version 2 of lZahnle et al.l (|2011 ). Since Zahnle 
et al.'s 100-bar lower boundary condition of a 1 x 10^^^ 
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Fig. 15. — Mole-fraction profiles for major species in a Ix solar 
composition model that assumes temperatures are constant at 1000 
K, JsTzz's are constant at 10^ cm^ s~^, and the planet is located 
at 0.1 AU from a Sun- like star, for a stellar zenith angle fixed at 
30°. These model results can be compared with isothermal models 
described in IZahnle eFall ll20Tll) . as well as with our HD 189733b 
and HD 209458b models results shown in Figs.|2j6] A color version 
of this figure is available in the online journal. 

mole fraction for specie s other than H, H2O , CO, CH4, 
NH3, N2, and H2S in the lZahnle et al.l ()201lD model leads 
to odd behavior for some of the minor constituents near 
the lower boundary, we simply assumed a constant con- 
centration gradient for minor species (including H and 
CO, which are minor species under these conditions at 
the 100-bar lower boundary), which means that these 
species will flow through the lower boundary at a max- 
imum possible rate. We also assumed equilibrium mole 
fractions for H2 and He at the lower boundary. 

A detailed comp arison of Fig. [15] with Fig. 4b of 
IZahnle et al.l (|2011l ) reveals some similarities and differ- 
ences between our models. The H2O mole fraction re- 
mains just below 10^^ in both models and is fairly un- 
perturbed by photochemistry. Atomic H is very abun- 
dant at high altitudes and greatly affects the chemistry of 
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molecular species in both models. Although our results 
for atomic H are qualitatively similar here, the quantita- 
tive results differ, likely because of differences in chem- 
ical reaction rate coefficients for hydrocarbon reactions. 
Methane has a sharp mole-fraction cutoff in both models 
at pressures of 10^"^ or 10^'' bar, although the location 
of the cutoff and shape of the CH4 profile is different 
in both models. This difference appears to be related to 
C2H4 chemistry, as the ethylene profiles are the most ob- 
vious difference between the models. The iZahnle et al.l 
([2011) results have C2H4 replacing CH4 at a relatively 
deep ~ 0.1 bar, causing a large "bite out" in the CH4 
profile from 0.1-0.001 bar before C2H2 takes over as the 
dominant hydrocarbon at high altitudes. The C2H4 in 
our model peaks at a somewhat higher altitude, and its 
abundance never approaches that of CH4, so methane is 
carried up to pressures of 10^** mbar before its eventual 
replacement by C2H2, HCN, and CO — note that both 
models agree that CO, HCN, and C2H2 will be the domi- 
nant carbon carriers at high altitudes. The differences in 
the hydrocarbon profiles and some portion of the H pro- 
files seem to stem from a different low-pressure limiting 
rate coefficient for H -|- C2H4 -I- M — >• C2H5 -I- M between 
our models. Due to a n inadvertent data-e ntry error, that 
rate coefficient in the IZahnle et al.l ()201lD model is three 
orders of magnitude smaller than the expression we have 
adopted, effectively cutting off an important C2H4 de- 
struction pathway. When this expr ession is changed t o 
better reflect literature values (e.g.. IBaulch et allll994D . 
the Z ahnle et al.l (|2011D results are more in line with ours 

(K. Zahnle, personal communication, 2010)^ 

Ano ther apparent difference between the IZahnle et al.l 
()2011|) models and ours is the CO and CO2 mole frac- 
tions at lower altitudes. Although the CO mole fraction 
remains at 10^^ to 10^"' at high altitudes in both mod- 
els, the drop off with decreasing altitude has a much 
different gradient in our model compared with theirs. 
Carbon monoxide is produced photochemically at high 
altitudes in both models, and CO then flows down to 
the lower atmosphere where in our model the large flux 
boundary con dition aflows it to be transported out of the 
system. Zahnle et al.l ()2011D do not discuss their domi- 
nant methane oxidation pathways in detail, so it is not 
clear how methane is converted to CO in their model, 
although they do mention that C -I- OH is important 
and they mention that they include three-body OH ad- 
dition reactions that can form C-0 bonds. In our model, 
methane is oxidized to CO through processes such as 

CH4 + H ^ CH3 + H2 
H + CH3 ^ ^CHa + H2 
H + 3CH2 ^ CH + H2 
CH + H2O ^ H2CO +H 
H2CO + H ^ HCO +H2 

HCO + H ^ CO + H2 

Net : CH4 + H2O + 4H ^ CO + 5H2; 

(15) 

however, this mechanism becomes ineffective once the 
H abundance drops significantly, and mechanisms like 
schemes (2) and (3) in section [3?2] take over but do not 
operate effectively at 1000 K. Therefore, in our model. 



diffusion con t rols th e CO profile at low altitudes. If the 
IZahnle et ail (|201lD model is fully converged, then it ap- 
pears that either the above CO production scheme (15) 
or alternative ones that produ ce the C-0 bond rem ain ef- 
fective at low altitudes in the IZahnle et al.l ()2011D model 
— perhaps due to their much larger low-altitude H abun- 
dance — allowing CO to be kinetically produced in the 
lower atmosphere. Whatever the cause, their increased 
CO abundance in the lower atmosphere allows an in- 
creased abundance of CO2 compared to our model, as 
CO and CO2 are kinetically linked. 

The qualitative results for the nitrogen-bearing species 
appear to be simil ar between our model and that of 
IZahnle et al.l (I20TI . Ammonia has a sharp cutoff at 
~ 10"'^ bar in both models, with HCN replacing NH3 as 
the dominant nitrogen carrier at higher altitudes. Molec- 
ular N2 is relatively unreactive in both models and has 
a constant mole fraction through out most of the atmo- 
sphere. The IZahnle et al.l ()2011D model apparently has 
an effective mechanism that converts HCN to N2 at high 
altitudes, whereas HCN is more efficiently recycled at 
high altitudes in our model. Similarl y, NO is produced at 
high altitudes more efficiently in the IZahnle et al.l ()2011[ ) 
model, or we have more efficient loss processes in our 
model. 

In general, our main conclusions from t his isothermal 
modeling exercise are similar to those of IZahnle et al.l 
()2011f ) in terms of the major carbon-, oxygen-, and 
nitrogen-bearing species and their relative importance in 
different regions of the atmosphere, with the exception of 
the C2H4 abundance, as discu ssed above. Howeve r, we 
disagree with the assertion of IZahnle et al.l (j2011l ) that 
such isothermal atmospheres are good analogs for extra- 
solar planets, and we feel the likelihood of soot formation 
on cooler exoplanets has been overemphasized as a result 
of Zahnle et al.'s reliance on isothermal models. Isother- 
mal model results do have some resemblance to those 
from models with more realistic temperature profiles 
(compare Figs. I3HS] with Fig. rT5|) . However, one critical 
difference is the presence of transport-induced quenching 
in the models with more realistic temperature profiles. 
Like the real HD 189733b atmosphere, the isothermal 
1000-K profile starts in the CH4 stability regime in the 
deep atmosphere but crosses over to the CO stability 
field in the upper atmosphere. However, unlike the real 
HD 189733b atmosphere, temperatures are low enough 
throughout the 1000-K isothermal atmosphere that inter- 
conversion between CO and CH4 (and between N2 and 
NH3, etc.) is inhibited everywhere. Therefore, whatever 
molecule is dominant at the base of the isothermal model 
will continue to be the dominant molecule throughout 
most of the atmosphere in the isothermal photochemi- 
cal models. Transport-induced quenching will not occur, 
and thus the abundance of species like CO will be un- 
derpredicted (as its only source will be photochemistry, 
not thermochemistry), and the CH4 abundance will be 
overpredicted as compared to kinetic-transport models 
that have higher temperatures at depth. 

This problem is particularly acute if the real planet has 
a temperature profile such that the composition would re- 
main in thermochemical equilibrium as the atmosphere 
transitions into the CO stability field, as seems to be 
the case for both HD 189733b and HD 209458b. In that 
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case, the CH4 — CO interconversion quench point wih 
occur while CO is dominant, and CO will be the main 
carrier of carbon in the stratospheres of these planets, 
in stark contrast to the isothermal model results. Even 
if the stratospheres become quite cold on these planets 
with more realistic temperature profiles, methane is not 
going to be as abundant as the isothermal models would 
indicate. Becaus e CH4 is the source of the complex hy- 
drocarbons in the lZahnle et al.l ()2011[ ) models (and in our 
own models) , a greater CO / CH4 ratio in the upper atmo- 
sphere implies fewer hydrocarbons and less chance of soot 
formation. The inability to form complex hydrocarbons 
is also exacerbated by the efficient hydrogenation of un- 
saturated hydrocarbons as the pressure increases, so that 
species like C2H2 and C2H4 are efficiently converted back 
to methane at pressures great er than 1 /xbar in our model 
(and probably would be in the lZahnle et al.l (|2011h model 
for more realistic considerations of the low-pressure lim- 
iting rate coefficient for H -I- C2H4 -I- M ->• C2H5 -I- M). 
Complex hydrocarbons are thus several orders of mag- 
nitude less im portant in our HP 1 89733b models than 
the isothermal iZahnle et al.l (|2011l ) models would sug- 
gest, and soot formation is much less of a factor. That is 
not to say that soot formation cannot still occur on these 
cooler exoplanets, as is discussed in section 13.41 partic- 
ularly if ion chemistry is conducive to the formation of 
heavy neutrals, but the soot precursors would be confined 
to high altitudes, and the overall column abundance of 
heavy hydrocarbons will be many orders of magnitude 
smaller than is suggested by the Zahnle et al. isothermal 
models. 

6. COMPARISONS WITH OBSERVED MOLE FRACTIONS 
ON HD 189733b 

Our model results can also be compared with the 
mole fractions inferred from the transit and eclipse ob- 
servations of HD 189733b (see Tabled]). The sugges- 
tion that C-, 0-, or N-bearing molecules are affect- 
ing the observational signature of HD 189733b during 
the primary transit and secondary eclipse came origi- 
nally from mid-infrared photometry from the Infrared 
Array Camera (IRAC) and from mid-infrared spectra 
from the Infrared Spectrograph (IRS ) on the Spitzer 
Space Telescove (iTinetti et all 120071: iGrillm air et al.l 



Beaulieu et al.|[2008l:lSing et al .1120011 Swain et al.l l2008bl: 



20081: iBeaulieu et al I l2008t iKnutson et al.l 120 07. 2008; 



Charbonneau et al .1 I2008t IPesert et al.l l2009r see also 



the th eore ti cal m odeling and mod el-data comparisons 
of'Fortncv & Mar lev 2007; Barman 20081 : iBurrows etHI 
[2OO81: .Madhusudhan fc Seager .2009). The mid-IR data 
have provided relatively loose constraints on molecu- 
lar abundances for HD 189733b (see Table [2]), whereas 
the near-infrared wavelength region is potentially more 
useful for constraining abundances because of the pres- 
ence of potentially diagnostic absorption features in 
this region. In fact, near-IR transit and eclipse 
data from the NICMOS instrument on the HST do 
seem to provide tight er constraints on molecular abun- 
dances (Swain et al. '2008b', '2 009al: iSing et al.l I200I 
iMadhus udhan & Scagcr 2009). However, the state 
of affairs in the analysis of transit and eclipse ob- 
servations is clearly not ideal, as contradictory data 
sets and/or analyses exist, both within the mid-IR 
and near-IR regions an d between the se two regions 
(e.g., IMadhusudhan fc Seager. ,2009; .Desert et al.l 120091: 



Gibson et al.ll2010l:lFortnev et al.ll2010[ ). If svst"ematic in- 
strument effects are not completely inhibiting the ability 
to acquire meaningful spectra l and photometric informa- 
tion fe.g.. iGibson et alJ 120101. we can at least attempt 
to compare our model results with the available observa- 
tions (see Tabled]) and discuss some general issues. 

Our 1 X -solar-composition kinetics-transport models 
for HD 189733b appear consistent with the loose con- 
straints supplied by the Spitzer IRS secondary eclipse 
observations of H2O and CH4 (Grillmair et al. 200l; 
IMadhusudhan fc Seager! 120091 ) . with the tighter con- 
straints from HST NICMOS secondary e clipse obser- 
vations of CO and CO2 as reported by I Swain et al.l 
(|2009al) . and with the CH4 and H2O abundances from 
the HST NICMOS transit o bser vations analyzed b y 
IMadhusudhan fc Seagerl f2009^ and lSwain et al.l (|2008b') . 
The models cannot, however, account for the much 
smaller CH4 mixing ratios inferred for the secondary 
ecli pse as opposed to the primary transit of HD 189733b 
(cf. iSwain et al.ll2008bl l2009at IMadhusudhan fc Seageii 
I2009f) . unless large temperature differences between the 
dayside and terminators can be maintained in the 1-100 
bar region and unless zonal transport does not homog- 
enize the CH4 abundance in the 0.1-1 bar regi o n. A s 
was first demonstrated bv iCooper fc ShowmanI (|2006f) . 
vertical transport-induced quenching is quite effective in 
the upper atmospheres of hot Jupiters, and CH4 is not 
readily converted to CO at dayside temperatures in our 
models. Some high-altitude CH4 is lost to HCN and 
C2H2 in our models, but there is still sufficient CH4 at 
transit-sensitive pressures that we would expect dayside 
observations to show CH4 mole fractions only a factor 
of a few smaller than at the terminators for the dayside- 
average and ter minator-average temperature profiles de- 
rived from the iShowman erall (2009) GCMs. If tem- 
peratures in the deep (1-10 bar) dayside atmosphere are 
much warmer than we have assum ed (e.g., compare the 
1-D 27r lFortnev et al.l ()2006all2010[ ) profile with the area- 
weighted dayside-average thermal profile shown in Fig. [2]) 
or if Kzz^s at the deep quench points are much smaller 
on the dayside than at the terminators, then this result 
could change. 

The fact that the predicted infrared eclipse depths 
from our dayside-average model seem to be uniformly 
too small compared to observations, combi ned with the 
low da y side CH4 abundance in f erred by ISwain et al.l 
(|2009al) : IMadhusudhan fc Seaged (pOOl . suggest a sce- 
nario in which the dayside- vs. -terminator temperature 
differences in the 1-100 bar regi on of HD 189733b are 
larger than is predicted from the IShowman et al.l (|2009( ) 
GCMs. Quenched CH4 abundances are m uch smaller in 
the model th at uses the warmer "27r" 1-D iFortnev et al.l 
I|2006all20ld[ ) thermal profile, for example (see Table p. 
The main problem with this scenario is that the radia- 
tive time constants at the CH4 quench level near 1-10 bar 
are expected to be greater than the zonal transport time 
scales, so that one would expect the temperatures to be 
homogenized longitudinally in this region. The strong 
zonal winds at would also be expected to help homoge- 
nize the CH4 abundance across longitudes, so that there 
would not be a big terminator vs. dayside difference, in 
conflict with the observations. 

As it stands, neutral C-H-O-N chemistry does not seem 
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TABLE 2 

Volume mixing ratios from HD 189733b observations and models 



Source 



H2O 



CO 



CH4 



CO2 



Secondary eclipse: Spitzer IRS 
Grillmair et al. (2008) 
Madhusudhan fc Seager (2009) 

Secondary eclipse: Spitzer IRAC 
^harbonneai^_et_aL ^2008) 
Madhusudhan fc Seager (2009) 

Secondary eclipse: HST NICMOS 
Swain et al. (2009a) 
Madhusudhan fc Seager (2009) 

Our dayside-average thermal-structure models at 0.1 



Ix solar 
lO"** to 0.1 

Ix solar 
(0.1-10) xlO"'' 

(0.1-1) xl0^4 
~ 1 X 10"" 



Ix solar 



(1-3) xlO-4 
(0.2-20) xlO"^ 



< 10" 



< 2 x 10" 



< 1 x 10 

< 6 xlO^ 



(7-700) xlO" 



(1-10) xlO-'' 
~ 7 X 10"* 



bar {i.e., for secondary-eclipse conditions) 



nominal 
-ftTzz = 10'^ < 
i^zz = lO^O 
nominal K, 



cm s 
lOx solar 



3.5 x 10 
3.4 X 10 
3.9 X 10 
3.4 X 10 



Our warmer "27r" ^brtne^_e^_ah 

nominal Kzz 

Primary transit: Spitzer IRAC 

Tinetti et al. (2007) 
Desert et al. (20^91 



'2006a . 2010) model 
3.4 X 10"'' 



4.5 X lO""' 

4.6 X 10-"' 
4.1 X 10"'' 

4.6 X 10-3 

at 0.1 bar {i.e.. 

4.7 X 10""' 



1.2 X 10 
2.9 X 10 
3.9 X 10 
6.1 X 10 



1.3 X 10 
1.3 X 10 
1.3 X 10 
1.3 X 10 



for secondary-eclipse conditions) 



3.7 X 10" 



9.7 X 10" 



~ 5 X 10""' 
upper limit 



~ 5 X 10-"' 
5 X 10-* to 0.1 
upper limit 



Primary transit: HST NICMOS 

Swain et al. (2008b) 
Madhusudhan fc Seager (2009) 
Sing et al. (2009) 

Our nominal K^z terminator-average model {i.e 

at 10-* bar 3.8 X 10-* 

at 10-3 bar 3.8 X IQ-* 

at 10-2 bar 3.8 X 10-* 

at 10-1 bar 3.8 X 10-* 



for 



1 X solar 
CO/H2O ~ 5-60 

~ 5 X 10-5 
lO--""' to 0.3 

primary-transit conditions) 

4.2 X 10-* 3.2 X 10-5 

4.2 X 10-* 4.3 X 10-5 

4.2 X 10-* 4.5 X 10-5 

4.2 X 10-* 4.6 X 10-5 



5.5 X 10 
4.5 X 10 
3.2 X 10 
2.0 X 10 



Note. — For the observations in which "Ix solar" or " ■ - ■ " (no data) is Ustcd, detailed constraints on 
the abundances were not presented. Our models assume Ix solar abundance, unless otherwise indicated. The 
"nominal" Kzz profiles are shown in Fig. [Tl 



to be able to account for a depletion of CH4 at I-IO"'^- 
bar pressures on the dayside of HD 189733b — where 
we expect transpor t-induced quench i ng to dominate — 
when we adopt the iShowman et all (2009) GCM-based 
temperature profiles. It is possible that sulfur photo- 
chemistry could contribut e to a loss of CH4 at these 
pressures on HD 189733b ()Zahnle et al.|[2009l[20TTI ). de- 
spite the low S/C ratio expected for solar-like compo- 
sitions, if catalytic destruction cycles are in operation. 
This process could only work effectively at low altitudes 
if the CH4 mole fraction were comparable to the H2S 
mole fraction (i.e., in situations where the quenched CH4 
abundance is small) because the solar ratio of S/C ^ 1 
will make it unlikely that products like CS and CS2 can 
make a dent in the methane abundance if CH4 is the 
dominant carbon carrier, and other potential products 
like C2Ha; hydrocarbons are unstabl e and will not r eplace 
methane at these pressures. In the IZahnle et afl (|20TTf ) 
isothermal models, for instance, CS and CS2 do not re- 
place CH4 at low altitudes because of the much greater 
CH4 abundance as compared with H2S. However, since 
CH4 is likely not the dominant carbon carrier on HD 
189733b and HD 209458b, sulfur chemistry might poten- 
tially play a l arger role in methane destruction than the 
IZahnle et al.l ()2011[ ) models indicate. Similarly, PH3 has 



weak bonds and could be an important source of atomic 
H at low altitudes that could affect CH4 chemistry, but to 
truly destroy CH4 in the lower stratosphere , some prod- 
uct like HCP must be able to compete with methane, 
which we think is unlikely given the thermodynamic sta- 
bility of HCP — sulfur photochemistry seems to have 
greater promise as a mechanism for destroying CH4 in 
the dayside lower stratosphere of HD 189733b. 

Although our predicted CO2 abundances for HD 
189733b are consistent with those from the HST 
NICMOS secondar y eclipse observations reported by 
ISwain et al.l (j2009a| ) , a much greater CO2 ab i indan ce is 
suggested from the IMadhusudhan fc Seageri (j2009D re- 
analysis of both these data and the Spitzer IRAC sec- 
ondary eclipse data. Carbon dioxide closely follows ther- 
mochemical equilibrium in our models. Some additional 
CO2 is produced at high altitudes from photochemistry 
in our models, but this amount is not sufficient to af- 
fect the column abundance at pressures greater than 
~10-^ bar. Therefore, if these larger observationally de- 
rived CO2 values are correct, then HD 189733b likely 
has an enhanced metallicity (see alsolFortnev et al.ll2010l : 
iLine et al.l 120101: IZahnle et al.ll20lH) . Enhanced metal- 
licities might also help explain the strong absorption in 
the 4.5 /Ltm band from Spitzer IRAC transit observations 
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HD 189733b terminator average, 1 0x solar 




Mole Fraction 



Mole Fraction 



Fig. 16. — Mole fraction profiles for several important species in dayside-average (left) and terminator-average (right) HD 189733b models 
that assume a lOx solar composition but otherwise have the same K^z profile, thermal structure, and background atmospheric grid as our 
nominal Ix solar HD 189733b models shown in Figs. [IHBl From comparisons between these figures, note the linear increase in the H2O 
and CO mole fractions, the quadratic increase in the CO2 mole fraction, the slight increase in the quenched NH3 mole fraction, and the 
slight decrease in the quenched CH4 mole fraction as the metallicity is increased from Ix to 10 X solar. A color version of this figure is 
available in the online journal. 



(IDesert et al.ll2()09(). A large CO/H2O ratio was invoked 
by iDesert et al.l ( 2009[ ) to explain this 4.5-/im feature, 
which would not be consistent wi th solar elemental ra - 
tios; however, as is described by iFortnev et al.l ()2010[ ). 
this feature could also be explained by the high CO2 
abundan ces that are predicted from enhanced- metallicity 
models (jLine et al.l l2010t iZahnle et afll2011[ ). A high 
metallicity for HD 189733b mi ght then explain both the 
iMadhusudhan fc Seage"rl ([2009') derived hig h CO2 abun- 
dances from the HST NICMOS data of iSwain et alj 
(I2q08bl l2009aD as we ll as the 4.5 ^m absorption seen 
by IDesert et'all (|2009D . 

Figure [16] illustrates how our HD 189733b model re- 
sults change if we assume a lOx solar metallicity, but 
otherwise keep the same Kzz profile and atmospheric 
grid (e.g., the T-P profiles and scale height stay the same 
despite the increased mean molecular weight and atmo- 
spheric opacity) as in our nominal 1 x solar model. The 
CO and H2O mole fractions increase by a factor of 10, 
but the CO2 mole fraction increases by a factor of '^100, 
and the CH4 mole fraction actually decreases when the 
metallicity is increased from Ix to lOx solar. The sensi- 
tivity of the CO2 ab undance to metallicity has been doc- 
umen ted previously (iLodders &: Feglevll2002t iLine et all 
[20Tnl : IZahnle et al l 120111 ) and is the natural result of 
C02's dependence on both CO and H2O. The decrease 
in the CH4 abundance with increased metallicity seems 
to result from the fact that the CH4 — )■ CO conver- 
sion schemes (2) and (3) in section 13.21 are more effec- 
tive when there is more H2O and CO present, such that 
CH4 follows its equilibrium curve to higher altitudes be- 
fore the quench point occurs, resulting in a decrease 
in the quenched CH4 abundance at higher metallicities. 
This higher-metallicity model seems at face value to be 
more consistent with the secondary-eclipse CO2 and CH4 
abundances inferred by IMadhusudhan k, Seagerl (|2009l ) 
and with the strong 4. 5-/im absorpt i on ob served in the 
transit observations of I Desert et al.l (|2009 ^. but the re- 
sulting water abundance may be higher than can be sup- 
ported by the HST NICMOS secondary eclipse data, un- 



less some additional absorption or scattering is muting 
the near-IR H2O absorption bands. 

In fact, the relatively low inferred H2O mixing ratio for 
HD 1 89733b from the HST NICMOS secondary eclipse 
data (ISwain et al.ll2009at IMadhusudhan fc Seagei1l2009( ) 
is somewhat disconcerting. Water is prevalent in both 
equilibrium and disequilibrium models for a wide range 
of metallicities. Only for greatly subsolar metallicities or 
greatly supersolar C/0 ratios can the H2O mole fraction 
be m uch less than 10"'^ in the atmosphere of H D 189733b 
(e.g. JLodders fc Fegleyll2Q0l ILine et ai1l2010l) . As is dis- 
cussed by numerous authors, water absorption is there- 
fore expected to dominate the infrared spectrum of this 
and most other transiting extrasolar giant planets. Given 
the photochemical stability of H2O and the likelihood 
of H2O being carried to as high altitudes as any other 
stable molecular constituent, it is unlikely that clouds 
or hazes could obscure H2O absorption features at near- 
infrared wavelengths without all other molecular features 
being obscured as well. However, where H2O does not 
strongly absorb, other molecules — including disequilib- 
rium molecules — might contribute opacity that could 
mask the expected signatures of H2O. As such, the 
observed apparent radius ratios between 1.66 /im and 
1.87 /xm in the HST NICMOS transit observations of 
iSing et all (2009) could indicate the presence of an addi- 
tional absorber rather than a lack of H2O. Alternatively, 
the contrast in the absorption bands could be washed out 
due to multiple Rayleigh or aerosol scattering, but such 
effects should mute the contrast for all near-IR absorp- 
tion features such that the relatively strong absorptions 
reported by Swain et al. (2008b, 2009a) could not be ex- 
plained. In any case, spectral models need to be updated 
to include the effects of disequilibrium constituents. 

These disequilibrium products will affect not only the 
observed spectra of HD 189733b, but the energy balance 
of the planet as well. A full discussion of the effect of 
disequilibrium constituents on the thermal structure of 
the planet are deferred to a future work. 
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7. COMPARISONS WITH OBSERVED MOLE FRACTIONS 
ON HD 209458b 

Although some infrared photometric and spectral 
data that can potentially help constrain the abun- 
dance of C-, 0-, and N-bearing molec ules on HP 
209458b have been acquired and analyzed (iBrown et al. 
2002; Richardson ct al. 2003, 2006, 200 7: IDeming et ah 
2p05a,b, 2007; Knutson ct al. 2008c jBeaulieu et al.ll200'a 
Madhusudhan & Scager 2009), actual reported con- 
straints on the molecular abundances are rare, due to 
signal-to-noise issues and to the fundamental degeneracy 
between temperatures and abundances that becomes a 
bigger problem for planets with stratospheric thermal 
inversions. Available constraints are listed in Table [31 
along with some results from our models. Our nominal 
HD 209458b dayside-average model, which is supposed 
to be representative of secondary-eclipse conditions, has 
a greater water mixing ratio, a smaller CH4 mixing ratio, 
an d a smaller CQ2 mix ing ratio than the values derived 
bv I Swain et al.l (j2009bD . Since the mid-IR Spitzer data 
themselves m ay not provide strict constraint s on the CH4 
mixing ratio ( 'Madhusudhan fc Seageii[2009D . this model- 
data discrepancy for CH4 may not represent a true prob- 
lem, unless the near-IR HST NICMOS data rather than 
the Spitzer IRAC and IRS data are the mai n source 
of the CH4 constraints from the ISwain et al.l (|2009bf ) 
analysis. In any case, it is hard to imagine how HD 
209458b could be very methane-rich, especially on the 
dayside, as the high stratospheric temperatures should 
drive the composition toward equilibrium, where CH4 
would then be a minor constituent. Larger eddy diffu- 
sion coefficients might help quench CH4 at greater abun- 
dances than we are seeing in our nominal model, but 
thermochemical reactions on the dayside would still act 
to limit the total column abundance that could be built 
up. Perhaps some other const i tuent i s masquerading 
as methane in the ISwain et al.l ()2009bf ) data. Hydro- 
gen cyanide is one possibility, although it, too, is pre- 
dicted to be less abundant on the warm dayside. Car- 
bon dioxide is another species that seems to be only 
loosely constraine d by the mid-IR photom etric data for 
HD 269459b ((Madhusudhan fc Seagei0OO9l). so our mod - 
els might not be out of line with the ISwain et"al] ()2009b[ ) 
data unless their CO2 constraints come largely from 
the HST NICMOS data rath er than the Spitzer pho- 
tometric data. However , both ISwain et al.l (12009b) and 
iMadhusudhan fc Seageii (|2009f ) constrain the H2O mix- 
ing ratio to be less than 10"'', which would seem to favor 
low metallicities for HD 209458b, and the constraints 
on other species in relati on to H2O seems to favor larg e 
(non-solar) C/0 ratios (Madhusudh an fc Seageii I2009D . 
We have yet to investigate these types of models. 

For the atmospheric composition at the terminators, 
even fe wer constraints exist for the molecular compo- 
sition. iBeaulieu et al.l (2009) do not provide any firm 
conclusions regarding molecular abundances and simply 
state that models that assume H2O mixing ratios of 10^"' 
to 10"'^ provide reasonable fits to the data, and that a CO 
mixing ratio of 10"'*, a CO2 mixing ratio of 10~^, and a 
CH4 mixin g ratio of 10~^ are allo wable. M odel-data com- 

firisons bv'Fortne v et al.l (|2010[ ) with the IBeaulieu et all 
009) Spi tzer IRAC data and the 24 /itm Spitzer MIPS 
point from lRichardson et al.l (|2006f) and H. Knutson (pri- 



vate communicatio n, 2009) do no t add a ny further con- 
straints, altho ugh iFortnev et al.l (1201 Of) note that the 
contrast in the IBeaulieu et al.l(|2009f ) mid-IR data is hard 
to reproduce with their equilibrium models. In particu- 
lar, absorption at 5.8 and 8 /^m is under-represented, or 
in general the contrast between the 3.6 /im and the 5.8 
and 8 /im apparent radii is under-represented by their 
models. As shown in section |4l and Fig. |9l the pre- 
dicted spectra from the disequilibrium models are simi- 
lar to spectra from equilibrium models, so disequilibrium 
chemistry of C, H, O, and N species cannot resolve these 
discrepancies. 

8. CONCLUSIONS 

Our kinetics-transport models for HD 189733b and 
HD 209458b demonstrate that disequilibrium processes 
such as photochemistry, chemical kinetics, and transport- 
induced quenching can dramatically affect the composi- 
tion of extrasolar-giant-planet atmospheres. Although 
thermochemical equilibrium can be maintained in the 
deep, hot tropospheres of these planets, rapid transport 
and the photochemistry that is initiated from the absorp- 
tion of ultraviolet photons from the host star can drive 
the composition away from equilibrium. The effects of 
these disequilibrium processes are more prominent for 
cooler exoplanets than for warmer exoplanets for two 
main reasons. First, the quench points for major species 
like CH4, CO, NH3, and N2 are dependent on tempera- 
ture: the colder the deep atmosphere, the less effective 
that thermochemical kinetics will be at maintaining equi- 
librium, and the deeper the pressure level at which the 
quenched species will depart from their equilibrium com- 
position. Since the major transport-quenched species 
CH4 and NH3 on HD 189733b and HD 209458b have 
equilibrium abundances that increase with decreasing al- 
titude, a deeper quench level corresponds to a greater 
mole fraction of the quenched constituent. Second, ther- 
mochemical kinetics will be less effective on planets with 
cooler stratospheres, such that disequilibrium composi- 
tions can be maintained. When stratospheric tempera- 
tures are high, such as is expected for HD 209458b and 
other planets with stratospheric thermal inversions — or 
that are very strongly irradiated in general — kinetic 
reactions are more effective at driving the composition 
back toward equilibrium, even in the presence of ultra- 
violet photons and transport-induced quenching. Dise- 
quilibrium processes are therefore expected to be more 
effective on the cooler HD 189733b than on the warmer 
HD 209458b. 

Although transport-induced quenching can modify the 
abundances of all species, the effects are typically more 
noticeable for the species that are not expected to be 
abundant in equilibrium at observable altitudes. When 
CH4-CO interconversion ceases to become kinetically ef- 
ficient, for example, both CH4 and CO will quench, but 
the quenching will be most obvious for CO on planets 
that have cooler upper atmospheres (e.g., brown dwarfs 
or giant planets far from their host stars) and most obvi- 
ous for CH4 on planets with warmer upper atmospheres 
(e.g., highly irradiated hot Jupiters). The location of 
the quench point — such as whether it is in the CH4- or 
CO-dominated regime — can also play a role determin- 
ing which molecule appears to be the most influenced by 
quenching. 
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TABLE 3 

Volume mixing ratios from HD 209458b observations and models 



Source H2O CO CH4 CO2 

Secondary eclipse: SpitzerlRAC (3.5, 4.6, 5.8, 8 fim), IRS (16 /im), and MIPS (24 /xm) 

Madhuaudhan fc Seager (2009'! < 10-** > 4 X 10"^ lO"** to 0.04 (0.02-70) X 10" 

Secondary eclipse: Spitzer IRAC, IRS and HST NICMOS 

Swain et al. ('2009b'l 8 X lO"'' to 1 X lO""* ••■ (0.2-2) X lO"* (0.1-1) X 10"' 

Our nominal Kzz dayside-average model at 0.1 bar (i.e., for secondary-eclipse conditions) 

nominal K^^ 3.4 X lO""* 4.7 X 10"'' 3.7 X 10"'' 5.7 X 10"* 

Primary transit: Spitzer IRAC 
jBeaulicu ct al. (2009) 1 X solar 

Our nominal western-terminator-average model (i.e., for primary-transit conditions) 

at 10-* bar 3.5 X lO"" 4.6 X lO"'' 6.3 X 10"^ 1.8 X 10"'' 

at 10-3 bar 3.5 X 10"* 4.6 x 10""' 7.3 x 10"^ 2.0 x lO"'^ 

at 10-2 bar 3.5 X lO"" 4.6 x 10""' 7.6 x 10"^ 1.9 x lO"'^ 

at 10-1 bar 3.5 X IQ-" 4.6 x 10"" 7.6 x 10"^ 1.2 X lO"'^ 



Note. — For the observations in which "Ix solar'' or " ■-■ " (no data) is listed, detailed eonstraints on 
the abundances were not presented. Our models assume Ix solar abundance, unless otherwise indicated. The 
"nominal" K^.^ profiles are shown in Fig. [Tl 



In our kinetics-transport models that include neutral 
carbon, nitrogen, and oxygen chemistry, we find that 
species like CO, H2O, CO2, and N2 are relatively un- 
affected by disequilibrium chemistry because of their 
strong bonds and/or efficient recycling. The vertical pro- 
files of these species are expected to follow equilibrium, 
except at very low pressures (P < 1 //bar), where photo- 
chemical processes can play a dominant role under some 
conditions. In contrast, the mole fractions of CII4 and 
NH3 are expected to be significantly enhanced {i.e., by 
several orders of magnitude) over equilibrium predictions 
in the bulk of the "photosphere" of these planets because 
of transport-induced quenching. In the upper strato- 
sphere, however, photochemical processes abruptly and 
effectively remove NH3 and CH4 on HD 189733b and HD 
209458b in favor o f photocheniical p roducts like HCN, 
C^Hp N, and C (cf. Lia ng et HIIMil IZahnle et al.l[2?M 
[201l[lLine et~l. 2010). 

As was first discuss ed by iLiang et"aLl ()2003[ ): lYelld 
(l200l :[G arcia Mufiozl ()2007l ). one major consequence of 
photochemistry on close-in transiting planets like HD 
209458b and HD 189733b is the production of huge quan- 
tities of atomic H due to catalytic cycles initiated by H2O 
photolysis. In our model, as with some models of thcrmo- 
spher ic photochemistry (e.g., lYellc 2004; Garcia Muno j 
l2007f ). H replaces H2 as the most- abundant constituent 
near the base of the thermosphere, and H will be the 
dominant neutral component throughout the bulk of the 
thermosphere of these planets. Our model differ s from 
those of Liang et al. (2003, 2004) and iLine et al.l (|2010l ) 
in the effectiveness of H2 destruction at high altitudes. 
Other significant photochemical products in our mod- 
els include HCN, unsaturated hydrocarbons like C2H2, 
atomic O, C, and N, and certain radicals like OH, CH3, 
and NH2. 

This disequilibrium chemistry will have some observa- 
tional consequences, as the disequilibrium products can 
affect the photometric and spectral signatures during 
the primary transit and secondary eclipse of transiting 
planets like HD 189733b and HD 209458b. Increased 



opacity from HCN, CH4, NH3, and under certain con- 
ditions C2H2 could be particularly important, especially 
for cooler planets like HD 189733b or for transit spectra 
that are influenced by the atmospheric transmission at 
the colder terminators. We recommend that investiga- 
tors who analyze transit and eclipse data should include 
the possible effects of nitrogen-bearing species like HCN 
and NH3 in their analyses. The additional disequilibrium 
species will also likely affect the thermal structure of the 
planets — a factor that was not considered in our present 
analysis — and future radiative models should include 
the effects of potential disequilibrium amounts of species 
like HCN, CH4, NH3, and C2H2. We note that spectro- 
scopic information on the hot bands of these molecules 
is needed to better predict the effects of these species 
on the observed spectra, and measurements of ultravio- 
let photoabsorption cross sections at high temperatures 
would improve the model predictions. 

We find that although we produce some benzene, com- 
plex nitriles, and other high-molecular-weight species 
in our HD 189733b and HD 209458b models, the pro- 
duction of these and other potential soot precursors is 
not very favorable with our currently adopted neutral 
carbon, oxygen, and nitrogen photochemistry, in part 
because CO rather that CH4 is the dominant carbon- 
bearing constituent at high altitudes in our models. In 
this respect, we disagree with the conclusi ons based on 
the cool, iso thermal model atmospheres of iZahnle et al.l 
([200I [20lTI ). On both HD 189733b and HD 209458b, 
we expect the CO O CH4 interconversion quench point 
to reside within the CO stability field, such that CO 
rather than CH4 is the dominant carbon-bearing con- 
stituent in the cooler upper atmospheres of these plan- 
ets — unless the eddy diffusion coefficient is extremely 
high (K,, > 10" cm2 s"i on HD 189733b). Since 
cool (< 1500 K) isothermal models are already quenched 
at the bottom boundary, where CH4 is dominant, they 
predict unrealistically high CH4 abunda nces on plan- 
ets like HD 189733b (jZahnle et al.l |2011[ ). The likeli- 
hood of soot formation on HD 189733b has therefore 
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been overemphasized in the iZahnle et alj ()201lD mod- 
els. On warmer planets like HD 209458b, complex hy- 
drocarbons will have particular difficulty surviving, as 
the high H abundance and efficient thermochemical ki- 
netics will greatly increase net destruction rates for these 
species. However, some soot formation may be possible 
on cooler planets like HD 189733b, particularly if one 
considers the possibility of ion chemistry that is initi- 
ated from extreme ultraviolet radiation or charged par- 
ticles from the parent star. On many planets with re- 
duced atmospheres within our own solar system, high- 
altitude hazes are produced from such processes (e.g. , 
West et al.lll98l l:IPrvor fc Hordl ll991t[Gera rd et a l.|[l995l: 
Wong et al.l(200 0, 2003; Fricdsoi^^ll2002; Wait e et al.l 
2007tllmanaka"fc Smithii2007tiVuitton et al.ii2007l . l2008[ ). 
We encourage the consideration of Titan-like ion chem- 
istry in future photochemical models of extrasolar giant 
planets. 

Our models of neutral carbon, nitrogen, and oxy- 
gen chemistry shed little light on the source of the 
absorber responsible for the stratospheric thermal in- 
versions inferred for HD 209458b and other exoplan- 
ets (e.g.,lBurrows ct al. 2007, 2008: Knutson et al. 2008, 
2009; Fortncv ct al. 2008a; Machalck et al. 2009 . 120101: 
Madhusudhan fc Seager .2009. .2010a : ,0'Donovan et al.l 
20101) . TiO was originally deemed to be a leading can- 
didate because it would be expected to condense on 
cooler, less- irradiated hot Jupiters while remaining in gas 
phase on warmer, more highly irra diate d hot Jupiters 
(iHubenv et al.ll2003l: iBurrows et all 120081 : iFortnev et al.l 
l2008al) : however. iKnutson et al.l ()2010[) argue that the 
lack of a convincing correlation between the incident 
flux and the presence of an inferred thermal inversion 
— as well as an observed apparent correlation between 
less chromospherically active stars and exoplanet thermal 
inversions — suggest that photochemistry might play 
a role, perhaps through destruction of the responsible 
absorber. In our models, we do not see any molecu- 
lar absorbers on HD 209458b that are not also present 
on HD 189733b. In fact, the cool stratosphere of HD 
189733b supports more disequilibrium molecular con- 
stituents than on HD 209458b. Atomic species are more 
prevalent on HD 209458b, and it is possible that some 
unidentified atom could be the culprit absorber, but that 
suggestion leads to a circular argument, as it is the high 
temperatures that cause atomic species to be more preva- 
lent in the first place. Chromospherically active stars 
have more fiares and produce more cosmic rays, and one 
might suggest that ion chemistry induced by charged par- 
ticles from the host star could produce a high-altitude 
haze that could obscure the culprit absorber — despite 
the fact that photochemical hazes tend to be absorp- 
tive th emselves and could thus heat the atmosphere (see 
iPont et al. 2008 for evidence of a Rayleigh-scattering at- 
mosphere at visible wavelengths on HD 189733b such 
that the Na absorption features are obscured, sug, gest- 
ing a high-altitude haze, and see lZahnle et al.ll2011l who 
champion soot formation at high altitudes). Alterna- 
tively, some sulfur p hotochemi cal product might con- 
tribute to absorption (jZahnle et al. 2009), or there could 
be a difference in dynamical mixing on different exoplan- 
ets. In any case, neutral C-N-0 chemistry seems unlikely 
to explain the difference in stratospheric thermal profiles 
on HD 189733b and HD 209458b because of the expected 



similar photochemical products on the two planets, but 
ion chemistry might contribute to the destruction of gas- 
phase absorber or to the formation of complex C- and 
N-bearing photochemical products that are ultimately 
condensible. 

The quantitative results of our models are very sensi- 
tive to factors such as the thermal structure, the eddy 
diffusion coefficient profile, and the planetary metallic- 
ity. Most of these parameters are not well constrained 
by observations. The eddy diffusion coefficients, in par- 
ticular, are difficult to obtain either theoretically or ob- 
servationally. If eddy diffusion coefficients in the lower 
stratosphere and upper troposp h ere w ere as low as are 
suggested bv lYoudin fc Mitchelll ()2010[ ). then transport- 
induced quenching of constituents would be much less 
important on HD 189733b and HD 209458b. We note 
that if the CH4 abundance and metallicity could be accu- 
rately constrained by observations, then one could place 
constraints on the value of Kzz at the CH4 — > CO quench 
point in the deep atmosphere. Similarly, one might sug- 
gest that the observed presence of heavy atomic species 
like O, C, and metals might provide some constraints on 
the eddy diffusion coefficients in the upper atmosphere, 
given that low K^z values of < 10* cm^ s~^ result in 
a very low homopause altitude in our models, making 
it difficult for heavy atomic species to be mixed into 
the thermosphere. Note, however, that this conclusion 
is based on hydrostatic models that do not include the 
hydrodynamic vertical wind, and it is likely that the hy- 
drodynamic wind could dominate in cases where Kzz is 
low, such that the concept of a homopause level has little 
meaning for close-in transiting exoplanets; heavy species 
could be dragged into the thermosphere regardless of Kzz 
in these situations, and our low Kzz models are likely not 
very realistic. 

Simple time-constant arguments suggest that advec- 
tion might effectively homogenize the longitudinal com- 
position such that the terminator-vs.-dayside differences 
in the abundance of quenched species like CH4 and NH3 
that are predicted by our simple 1-D models would be 
washed out at pressures greater than a few mbar. There- 
fore, observations that indicate a difference in the col- 
umn abundance of these species during transit as com- 
pared with eclipse observations could provide clues to the 
rate of horizontal transport in the lower stratosphere. In 
the upper stratospheres at pressures less than < 1 mbar, 
photochemical time constants likely dominate, suggest- 
ing that the composition could be highly variable with 
longitude at high altitudes. One important outcome is 
the potential conversion of HCN and C2H2 back to CH4 
at high altitudes on the nightside hemisphere. When that 
high-altitude methane flows back around to the dayside, 
it may strongly absorb stellar radiation, with some inter- 
esting radiative consequences, including potential emis- 
sion in the band of methane that has some strong 
lines that could potential absorb at high alti tudes if CH4 
is ab undant enough (see the observations of iSwain et al.l 
l2010f ). The quantitative details need to be worked out to 
determine if this process is viable. 

The qualitative details of the photochemistry described 
here will be similar to other hot Jupiters and hot Nep- 
tunes. On very cool planets like GJ 436b, CH4 will dom- 
inate over CO, with a resulting increase in the produc- 
tion of photochemical products like hydrocarbons and 
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nitriles. However, CH4 will not be efficiently removed 
from the lower atmosphere of GJ 436b with our stan- 
dard neutral C, O, N chemistry, and we cannot ex- 
plain the low inferred CH4 abundance on this planet 
(jStevenson et a l. 2010). Sulfur photochemistry could po- 
tentially help (.Zahnlc et al. ,2009, 201 1.). but only if the 
total sulfur mole fraction approaches the CH4 mole frac- 
tion (i.e., for S/C ratios ^ solar) or if as-yet- unidentified 
catalytic cycles can permanently convert the methane 
into other hydrocarbons or C-0 bearing species. 

We eagerly await future results from the upcoming 
James Webb Space Telescope (JWST)^ whose broad 
wavelength coverage and other spectral capabilities 



should help reduce current uncertainties in the determi- 
nation of atmospheric composition on extrasolar giant 
planets. 
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